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ABSTRACT 

We present samples of starburst galaxies that represent the extremes discov- 
ered with infrared and ultraviolet observations, including 25 Markarian galaxies, 
23 ultraviolet luminous galaxies discovered with GALEX, and the 50 starburst 
galaxies having the largest infrared/ultraviolet ratios. These sources have z < 
0.5 and cover a luminosity range of ~ 10 4 . Comparisons between infrared lu- 
minosities determined with the 7.7 /im PAH feature and ultraviolet luminosi- 
ties from the stellar continuum at 153 nm are used to determine obscuration in 
starbursts and dependence of this obscuration on infrared or ultraviolet lumi- 
nosity. A strong selection effect arises for the ultraviolet-selected samples: the 
brightest sources appear bright because they have the least obscuration. Obscu- 
ration correction for the ultraviolet-selected Markarian+GALEX sample has the 
form log[UV(intrinsic)/UV(observed)] = 0.07(±0.04)M(UV)+2.09±0.69 but for 
the full infrared-selected Spitzer sample is log[UV(intrinsic)/UV(observed)] = 
0.17(±0.02)M(UV)+4.55(±0.4). The relation of total bolometric luminosity L ir 
to M(UV) is also determined for infrared-selected and ultraviolet-selected sam- 
ples. For ultraviolet-selected galaxies, log L ir = -(0.33±0.04)M(UV)+4.52±0.69. 
For the full infrared-selected sample, log L ir = -(0.23±0.02)M(UV)+6.99±0.41, 
all for Li r in L and M(UV) the AB magnitude at rest frame 153 nm. These 
results imply that obscuration corrections by factors of two to three determined 
from reddening of the ultraviolet continuum for Lyman Break Galaxies with z > 
2 are insufficient, and should be at least a factor of 10 for M(UV) ~ -17, with 
decreasing correction for more luminous sources. 
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Introduction 



The formation and evolution of gal a xies in the Univer se is now being observed to red- 



shifts z ~ t 
wavelengths 



e.g. 



Bouwens et al.l 12009c lOuchi et al.l 120091 ) using photometric imaging to 



1 //m in search of "Lyman Break Galaxies" (LBGs; Steidel et al. 1999). 
These observations reveal the rest frame ultraviolet luminosity of star-forming galaxies and 
allow the measurem ent of changes in th e star formation rate density of the Universe as a 
function of redshift (IMadau et al.lll998l ). (Because such star-forming galaxies are observa- 
tionally characterised by indicators of young, massive, short-lived stars, we refer to them in 
the remainder of this paper as "starburst" galaxies.) 

Observations of LBGs indicate an increase in the characteristic luminosity of starburst 
galaxies for < z < 2, and then a decr ease in luminosity for z > 2 until the observational 
limit at z ~ 8 ( IReddy and Steidell 120091 ) . Is this evolution of starburst galaxies determined 
correctly based on LBGs? The most uncertain factor in quantifying the luminosities is 
determining an accurate correction for obscuration by dust, because this obscuration can be 
extreme in the rest frame ultraviolet. 

The importance of dust in the early universe is emphasized by the result that the most 
luminous galaxies known have most of their luminosity emerging as reradiation by dust of 
intrinsic luminosity initially generated at much shorter wavelengths and absorbed by the 
dust. Initially termed " Ultraluminous Infrared Galaxies" (ULIRGS; Soifer et al. 1987) when 
discovered with the Infrared Astronomical Satellite (IRAS), such galaxies ar e found in large 



numbers in surveys at mid-infrared wavelengths at all redshifts to z ~ 3 (e.g.lLe Floc'h et al. 



2005 



Dev et al.l 120081) and millimeter wavelengths (jChapman et al. 1120051 ; iPope et al 



2008 



Menendez-Delmestrel boo"^. The evolu tion of starbursts among ULIRGs has been tracked 
to z ~ 3 ( jWeedman and Houck 2008 ) with the Spitzer Space Telescope using redshifts 



measured with The Infrared Spectrograph on Spitzer (IRS; Houck et al. 2004). Because of 
the inability to measure redshifts of dusty starbursts with z > 3 using the IRS, these have 
not yet been found at redshifts as high as LBGs. 

Combining our understanding of dusty galaxies dominated by infrared luminosity with 
galaxies discovered by ultraviolet observations is essential to a full understanding of the 
formation and evolution of galaxies (as well as the formation and evolution of dust). The 
availability of IRS spectroscopy for the same galaxies which are detected in the rest frame 
ultraviolet by the Galaxy Evolution Explorer (GALEX, Martin et al. 2005) allows a direct 
comparison between the infrared luminosity from dust and the ultraviolet luminosity from 
hot stars. 



Our initial comparison in this way for 287 starburst galaxies having both Spitzer spectra 
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and GALEX ultraviolet measures found large differences in the sta r formation rates (SFRs) 



dedu ced from infrared compared to ultraviolet measures of SFR (jSargsyan and Weedman 



20091 ). These differences are attributed to extreme obscuration of the ultraviolet luminosity. 



Our previous sample arose primarily from selection in the infrared, chosen as starbursts 
having IRS spectra available. In the present paper, we use IRS spectra for samples of 
ultraviolet-selected starburst galaxies to study obscuration effects in sources with similar 
selection criteria as LBGs, in order to understand how conclusions depend on selection 
criteria. Our objective is to use observations in the infrared and ultraviolet to make an 
indep endent test of dust correc tions for starbursts detected as LBGs at high redshift, such 
as in iReddy and Steidell ( 120091 ). and to determine reliable corrections for dust obscuration 
which are independent of selection technique. We make measures only of the obscuration in 
order to avoid uncertain assumptions regarding transformation of luminosities to SFRs. 

In section 2, below, we describe the sample selection which includes 25 Markarian galax- 
ies, 23 ultraviolet luminous galaxies discovered with GALEX, and the 50 starburst galaxies 
having the largest infrared/ultraviolet ratios. In section 3, we present the infrared spectra of 
all 98 sources and compare the characteristics of infrared and ultraviolet spectra in section 
4. Results are used in section 5 to develop an empirical method for determining obscuration 
corrections as a function of M(UV), and in section 6 to compare bolometric luminosities Lj r 
with M(UV). 

Our most important result (section 7) indicates that the bolometric luminosities deduced 
for starburst galaxies using LBGs at z > 2 have been underestimated, because obscuration 
corrections have been significantly underestimated for starbursts. This means that SFRs 
have also been underestimated at high redshifts, and we propose a new observational test to 
verify these conclusions. 



2. Selection of Starburst Galaxies 

For the various comparisons which follow in this paper, we define three independent 
samples of starburst galaxies which have been selected using either ultraviolet or infrared 
criteria. These are the "Markarian" sample, the "GALEX" sample, and the Spitzer sample. 
These samples are then compared using various spectroscopic parameters and conclusions are 
reached regarding how measures of dust obscuration and SFRs depend on sample selection 
criteria and source luminosities. 
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2.1. Markarian Starburst Galaxies 



The first sur vey for galaxies selected on the basis of ultraviolet brightness was the 
Markarian survey (jMarkarianlll967l ) . which detected "galaxies with an ultraviolet continuum" 
based on their apparent brightness in objective prism spectra extending to wavelengths > 340 
nm. Although the Markarian survey was not homogeneous for flux limits or selection criteria, 
it produced large samples of galaxies which subsequently proved important for classification 
of ultraviolet-luminous galaxies. For examp le, the prototype "star burst galactic nucleus" was 
defined by Markarian 538, also NGC 7714 ( jWeedman et al.lll98ll ). The first set of starburs t 
galaxies arising from this definition was the list of 50 Markarian galaxies in iBalzanol ( 119831 ). 
selected using optical spectroscopy for measurement of emission line widths to distinguish 
starbursts from active galactic nuclei (AGN). 

Because the selection of these Markarian starburst galaxies was based only on optical 
criteria, these galaxies provide an important comparison sample to starbursts selected in 
the infrared (section 2.3). (It was subsequently found that most Markarian galaxies would 
also be detected by the Infrared Astronomical Satellite (IRAS), but IRAS had not flown 
when the Markarian lists were compiled.) For this reason, we present in this paper new IRS 
observations from our Spitzer program 50834 of 25 Markarian galaxies selected from the 
Balzano (1983) list, and also include archival observations of NGC 7714. 



2.2. GALEX Starburst Galaxies 



The availability of the GALEX all sky survey enables the identification of nearby galaxies 
(z < 0.3) which are similar in ultraviolet luminosity to the dist ant LBGs; these local galaxies 
have been called ultraviolet luminous galaxies, or UVLGs (IHeckman et al.l 120051 ) . They 
represent the extremes in luminosity of ultraviolet-selected local starbursts, and they are 
sufficiently bright that various other data are available. 

A number of these UVLGs selected with GALEX have IRS observations available 
(Spitzer archival program 40444, L. Armus, P.I.). We have made new extractions of the 
IRS spectra for these sources and refer to them as the GALEX sample. 



2.3. Highly Obscured Spitzer Starburst Galaxies 



Our previous paper (jSargsyan and Weedmanll2009i ) compared IRS spectra with GALEX 
data for 287 starburst galaxies with z < 0.5, sources chosen only by having both IRS spec- 
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tra and GALEX fluxes. In that paper, the SFR measured from the luminosity of the 7.7 
fim polycyclic aromatic hydrocarbon (PAH) feature, SFR(PAH), was compared to the SFR 
measured from the rest frame ultraviolet luminosity at 153 nm, SFR(UV). The relation 
found was log [SFR(PAH)/SFR(UV)] = 0.531og [i/L„(7.7 /mi)] - 21.5 for SFR in M yr^and 
jUm) in erg s _1 . This relation indicates a much higher SFR as measured in the in- 
frared, which indicates very severe extinction in the ultraviolet for these starbursts, implying 
the escape of < 1% of the ultraviolet luminosity for the most luminous starbursts. 

To compare with the Markarian and GALEX ultraviolet-selected starbursts as one ex- 
treme of selection, we desire a comparably sized sample representing the other extreme of 
infrared selection. For this sam ple, we utilize the infrared-selected starbursts tabulated in 
Sargsyan and Weedmanl (120091 ) discovered in Spitzer flux-limited surveys and having the 
most extreme values of SFR(PAH)/SFR(UV), those with log [SFR(PAH)/SFR(UV)] > 2. 
This gives a sample of 50 sources, which we call the Spitzer sa mple. Original IRS observa- 
tions of this Spitzer sample are in IWeedman and Houckl ( 120091 ) . taking sources from a flux 
limited survey to f„(24/tm) > 10 mJy with the Multiband Imaging Photometer for Spitzer 
(MIPS, Rieke et al. 2004), and in archival program 40539 (G. Helou, P.I.), which derives 
sources from various flux limited surveys with MIPS to 5 mJy. 

This Spitzer sample also provides a low redshift sample that approaches in infrared 
luminosity the high luminosity, optically faint, and yery dusty starbursts discov ered with 
Spitzer at z ~ 2 (source list a nd original references in I Weedman and Houck plus new 



sources m 



Desai et al.l (120091 )). which have been defined as "Dust Obscured Galaxies", or 



DOGs (IDey et al.ll2008l ). In the same fashion, therefore, that our ultraviolet selected samples 
represent local LBGs, the Spitzer sample represents local DOGs. (Luminosity comparisons 
between DOGs and the local sample are shown in Figure 16.) 



3. Observations and Results 

3.1. Spitzer Infrared Spectra 

The infrared spectra used for our analysis were all obtained with the IRS using the 
Short Low module in orders 1 and 2 (SL1 and SL2) and the Long Low module in orders 1 and 
2 (LL1 and LL2). Combining these orders gives total low resolution spectral coverage from 
~5 /mi to ~35 /mi. All spectra which are used in this paper are presented below, although 



1 The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded 
by NASA through the Jet Propulsion Laboratory and the Ames Research Center. 
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spectra of 7 sources were originally given in IWeedman and Houckl (120091 ) . 



An IRS observation in one order consists of two independent spectra obtained at two 
different nod positions on the slit. Each spectrum is derived from several different images 
taken at the same nod position. All images with the source in one of the two nod positions 
were coadded to obtain the image of the source spectrum. The background was determined 
from coadded background images that added both nod positions with the source in the other 
slit (i.e., both nods on the LL1 slit provide a background image when the source is in the 
LL2 slit). The difference between coadded source images minus coadded background images 
was used for the spectral extraction, giving two independent extractions of the spectrum for 
each order. These independent spectra were compared to reject any highly outlying pixels 
in either spectrum, and a final mean spectrum was produced. 

The extraction of one dimensio nal spectra from the two dimensional imag es was done 
with the SMART analysis package (IHigdon et al.l 12004 ; iLebouteiller et al.ll2009l ). beginning 
with the Basic Calibrated Data products, version 18, of the Spitzer flux calibration pipeline. 
Sources in the Spitzer sample and GALE X sample are faint and a lways unresolved, so that 
use of the SMART Optimal Extraction in ILebouteiller et al.l ( 120091 ) improves signal to noise 
(S/N). The brighter Markarian sources, sometimes slightly resolved, were extracted with the 
"tapered column" extraction which uses an extraction window that is wider in the direction 
perpendicular to dispersion. 

That some Markarian starburst galaxies are slightly resolved at the scale of the IRS slits 
is noticed when the overlapping portions of the SL1 and LL2 spectra show a larger flux in the 
wider LL2 slit (the SL1 slit is 3.7 "wide, and the LL2 slit is 10.5 "wide). When the SL1 and 
LL2 flux densities are different in the overlapping wavelengths, the spectra are matched, or 
"stitched", by raising the SL1 by the factor needed to match LL2. These "stitching" factors 
are given in Table 1. 

Figure 1 shows spectra of all starbursts in the Markarian sample, Figure 2 all starbursts 
in the GALEX sample, and Figure 3 all starbursts in the Spitzer sample. Twelve AGN 
found among the Markarian and GALEX samples (section 4.2) are shown in Figure 4. In 
all Figures, spectra are truncated at 25 ixm in the rest frame because of the absence of any 
significant features in the low resolution spectra beyond that wavelength. Illustrated spectra 
are not smoothed. For convenience of illustration, spectra are arranged by spectral slope; 
i.e by the ratio i v {2& //m)/f„(7.7 jum). Spectra in Figures are normalized to 1 mJy at 7.7 
jim and offset for illustration. Observed spectroscopic data for all sources are given in Tables 
1 and 2 (Markarian sample), 3 and 4 (GALEX sample), and 5 and 6 {Spitzer sample). 



Average spectra for all starbursts in the 3 samples are shown in Figure 5. These averages 
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were produced by averaging the individual normalized spectra of starburst galaxies in Figures 
1, 2, and 3. After averaging, the average spectra were boxcar-smoothed to 0.2 /iin for 
illustration in Figure 5. 



3.2. GALEX Ultraviolet and SDSS Optical Data 



GALEX ( Martin et al.l I2OO0] ) obtains images i n the ultraviolet at wavelength s 134-179 



nm (FUV) and 177-283 nm (NUvjE As discussed in lSargsvan and Weedman] J2009h . GALEX 



ultraviolet-selected starbursts and Spitzer infrared-selected starbursts generally agree well in 
position, to accuracies much better than the observing apertures, which are of similar sizes. 
Point source GALEX images have full width at half maximum (FWHM) of 4.2" for FUV and 
5.3" for NUV (Morrissey et al. 2005). For the Spitzer IRS, spatial resolution is wavelength 
dependent, so the FWHM of an unresolved source at the mid-infrared wavelengths used for 
our spectroscopic analysis is 3"< FWHM < 6". 

The observed flux density at a rest frame FUV wavelength of 153 nm is determined by 
using the GALEX observed flux densities at FUV (effective wavelength 153 nm) and NUV 
(effective wavelength 227 nm) to fit a power law continuum between these wavelengths, and 
then interpolate to the flux density at observed wavelength (l+z)153 nm using the resulting 
power law index. This f iy [(l+z)153 nm] is the value of f„(FUV) listed in Tables 1, 3, and 5. 
When we refer in this paper to ultraviolet luminosities or absolute magnitudes M(UV), we 
mean M(153 nm). 

For comparison with the IRS spectra, we also use optical spectra from the Sloan Digital 
Sky Survey (SDSS, Gunn et al. 1998). We classified all such spectra to determine if the 
Balmer lines show broad wings characteristic of type 1 AGN, or have unusually strong 
[OIII]/H/3 ratios, which may indicate type 2 AGN. Sources with H/3 comparable in strength to 
[OIII] A 4959 without broad H/3 wings are classified as starbursts. These SDSS classifications 
result in some of the AGN discussed in section 4.2. 



2 GALEX data used in this paper were those available in August, 2009, at the MultiMission Archive at 
Space Telescope Science Institute |http: / / archive. stsci.edu/| 
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4. Infrared Spectroscopic Characteristics of Starbursts 



4.1. Infrared Luminosities and PAH Strengths 

Starb urst galaxies show very similar spectra in the mid-infrared, characterised by stron g 



PAH emission features (e.g. iRigopoulou et al.ll2000l ; iGenzel et al.lll998c iBrandl et al.ll2006l ). 
The exceptional similarity among such spectra is well illustrated in Figures 1, 2, and 3. 
Starburst spectra can be separated from spectra containing a contribution to dust emission 
from an AGN by using the equivalent widths (EW) of the PAH features. The best PAH 
feature for such quantitative use in starburst spectra observed with the IRS over all redshift 
ranges < z < 3 is the 6.2 /im feature, because this feature is spectrally isolated and does 
not require deconvolution from adjacent features. 

Because of the strength of the PAH features, luminosities of these fe atures provide a mea- 
sure of starburst luminos ity which relates to the star formation rate (e.g. iForster Schreiber et al 
2004 ; IBrandl et al.l 120061 ) . The PAH luminosity parameter which we have defined and used 
previously is vL u (7 .7 (j,m) , where L v (7 .7 urn) is deter mined from the flux density f u (7.7 //m) 
at the peak of the 7.7/im feature (IHouck et al.l 120071 ). This parameter was chosen to enable 
comparisons with starbursts at redshifts z ~ 2, for which the 7.7 /im feature is the strongest 
spectral feat ure and whose peak flux den sity is the simplest parameter to measure in sources 
of poor S/N ( jWeedman and Houckj|2008l ). Measurement of a peak f,(7.7 /im) also avoids the 
uncertainties of deconvolution and definition of the underlying continuum for the broad PAH 
features, which can lead to different measures of PAH flux even from the sa me spectrum. 



Altho ugh such deconvolutions are valuable for sources of sufficiently high S/N ( jSmith et al. 



20071 ). they are difficult to implement in sources at high redshift with poor S/N and few 
visible PAH features. 

An EW(6.2 /im) > 0.4 /im has been found empirically to be a good discriminant for 
"pure" starbursts; sources with a PAH feature of this strength invariably have optical spectra 
classified as starbursts. We examine this criterion for the present samples in Figure 6, where 
EW(6.2 /im) is compared with source luminosity determined from z/Lv(7.7/im). In this 
Figure, the starbursts in Figures 1, 2, and 3 are distinguished from the AGN in Figure 
4 (classification of the AGN sources is discussed below in section 4.2). The 3 GALEX 
sources with starburst classification but EW(6.2 /im) < 0.4 are especially interesting and 
are discussed in section 4.2 regarding whether they are really starbursts or have an AGN 
contribution. 

Our new results provide an independent confirmation of the EW criterion adopted 
previously to define pure starbursts. It is expected that the Spitzer sample should have 
EW(6.2 /im) > 0.4 /im , because this was a criterion for inclusion in the tabulation of 
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Sargsyan and Weedmanl (120091 ) . from which we drew this sample. It is notable, however, 



that the Markarian starbursts all satisfy the criterion of EW(6.2 iim) > 0.4 /im . This is in- 



dependent confirmation of the EW criterion, because the lBalzanol ( 119831 ) sample of Markarian 
starbursts was defined using only optical criteria. The four Markarian sources with EW(6.2 
jum) < 0.4 /im now appear with SDSS spectra to be AGN, so the original classification was 
incorrect, overlooking weak, broad wings on the hydrogen lines. 

We note, however, that two GALEX sources with a starburst classification from both 
SDSS and IRS spectra have EW(6.2 /im ) < 0.4 /im . These are sources 3 and 16 in Figure 
2B (the low EW measurement for source 13 is not reliable because the feature is so close 
to the end of the spectrum). In previous studies, such weak PAH strengths in sources 
without AGN were observed only in blue compact dwarfs (BCDs; Wu et al. 2006). These 
two sources have luminosities greatly exceeding BCDs, however, so they are potentially an 
important class of luminous starburst with weak PAH features. As discussed in efforts to 
understand BCD spectra, such sources could arise if the ionization radiation is unusually 
hard, such as in a very young starburst, or if the HII region is density bounded and not 
surrounded by molecular clouds, or in circumstances of low metallicity. 

No such sources with low EW(6.2 /im) are encountered among the Markarian sample, 
and the presence of several other AGN within the GALEX sample (section 4.2) raises the 
question of whether these objects also have a strong dust continuum and weak PAH because 
of the presence of an unrecognized AGN. Because of this possibility, we can at present only 
call attention to the unusual nature of these sources but cannot confirm their nature. 

Figure 6 also shows the luminosity range of our samples, a factor of 10 4 , and illustrates 
that all three samples overlap in luminosity. Especially notable is that the ultraviolet- 
selected starbursts provide a continuous range of luminosities, typically lower luminosities 
from the Markarian sample and higher from the GALEX sample, but they overlap in both 
EW and luminosity parameters. That these samples also overlap in z^L !/ (7.7//m) with the 
Spitzer sample of highly obscured starbursts is evidence that these very different selection 
techniques have nevertheless returned samples of similarly luminous starbursts and with 
similar spectroscopic nature regarding strength of the PAH star formation indicators. 



4.2. Sources Classified as AGN 

When comparing ultraviolet and infrared luminosities in order to study starbursts, it 
is important that AGN are not included in the samples. An AGN can raise the ultraviolet 
luminosity and also raise the level of infrared dust continuum. We are not considering details 
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of AGN in the present paper, but after definition of the samples described in section 2, it 
was found that spectra showed AGN in a few cases. Because we have IRS spectra of those 
sources, we identify them and show the spectra, but we do not use them subsequently in our 
analysis of starbursts. 

In Tables 1 and 3, sources in the Markarian and GALEX samples which are classified by 
us as AGN are noted. These spectra are shown in Figure 4, and these sources are not included 
as starbursts in the various Figures. For the five Markarian AGN (Mkn 93, 149, 193, 764, 
and 863), the classification as AGN arises from SDSS spectra which sho w weak, broad wings 



to the Balmer lines. These were overlooked in the initial classification of lBalzanol ( 119831 ) from 
which we drew the sample. All IRS spectra for these AGN except for Mkn 193 (source 16 in 
Table 1) show strong PAH features. This means that the infrared spectrum is characterized 
by a circumnuclear starburst, and the PAH features could be used to measure the luminosity 
of this starburst. The ultraviolet luminosity may be contaminated by the AGN, however, so 
we do not use these sources to compare ultraviolet and infrared luminosities. Mkn 193 also 
shows strong [SIV] 10.5 /mi emission, which is an indicator of the hard ionizing radiation 
from an AGN. 

For the GALEX sources in Table 3, we classify 7 as AGN. This classification includes any 
source which could be an AGN based on either optical or infrared spectra. Such conservative 
criteria are designed to assure that we do not include sources with ultraviolet luminosity 
arising in part from an AGN. In most cases, the AGN classification arises because of weak 
PAH 6.2 /mi together with absorption in the silicate feature at ~ 10 /mi . This combination 
is interpreted as evidence of an AGN with hot dust close to the AGN whose emission raises 
the level of continuum beneath the PAH and is also absorbed in the silicate feature by 
intervening cooler dust (Imanishi 2007). 

Five GALEX AGN are classified in this way with IRS spectra (sources 1, 2, 4, 10, and 
23 in Table 3 and Figure 4). There is some ambiguity in this classification, however. Sources 
3 and 16 in Table 3 are included in Figure 2 as starbursts despite having weak EW(6.2 /mi), 
because they show no evidence of silicate absorption. Sources 3 and 16 have high [OIII]/H/? 
ratios from SDSS spectra but no evidence that the lines are broad as in type 2 AGN. Such 
high ratios also arise in BCDs, for example, because of harder ionizing radiation, so we do 
not exclude these as starbursts. The question of emission line ratios and ionizing radiation is 
considered more in section 4.3, below. GALEX sources 12 and 15 from Table 3 are classified 
as AGN both because of SDSS spectra and because of strong [SIV] 10.5 /mi in the IRS 
spectra. 

These Markarian and GALEX AGN are not included in the remaining discussion, be- 
cause we do not want to confuse conclusions regarding starbursts by including sources whose 
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spectroscopic parameters may be influenced by an AGN. 



4.3. Dust and Ionizing Radiation in Starbursts 

Mid-infrared spectra of starbursts such as those in Figure 5 contain data for several 
independent indicators for the luminosity of the starburst. These are the luminosity of the 
PAH emission features discussed above, the luminosity of the dust continuum underlying 
the PAH features, and the luminosity of the emission lines. These three indicators measure 
different characteristics of the starburst. The PAH emission is excited by photons of various 
energies penetrating the photodissociation r egion (PDR) at the boundary between the HII 



region and the surrounding molecular cloud (IPeeters et al.l 120041 ) . Such photons are not the 



most energetic photons which have ionized the HII region, because such energetic photons 
would destroy the PAH molecules. The dust continuum arises from dust intermixed in the 
HII region and heated primarily by the hot stars, and the emission lines arise within the HII 
region from ionizing photons arising in the hot, young stars. 

Differences among these three parameters (PAH, dust, emission lines) give information 
about differences in the nature of different starbursts. Some notable differences are seen 
when comparing the three samples of starbursts considered herein. These differences are first 
illustrated in Figure 5, which shows average spectra for each sample. The well established 
dominance of PAH features in starburst spectra is again confirmed in the three samples, as 
shown for individual sources in Figures 1, 2, and 3. There are no measurable differences 
among relative PAH strengths for the three samples. There are, however, differences in the 
relative strengths of the emission lines and the PAH features. 

The relative strengths of the emission lines compared to PAH gives a measure of relative 
luminosities from the HII region compared to the surrounding PDR where the PAH features 
arise. Relative intensities of the emission lines themselves provide a diagnostic for the hard- 
ness of ionizing radiation in the HII region. The most important emission lines observed in 
the low resolution IRS spectra are labeled in Figure 5. 

The [Nell] 12.8 /im emission line is blended with an adjacent PAH feature so the actual 
measurement of its flux is uncertain in these low resolution spectra, but the shape and 
relative flux of this blended feature does not change among the three samples. The total flux 
of the blended feature is dominated by [Nell], so we can conclude that the HII luminosity 
relative to the PDR luminosity is similar among the samples. Compared to the fluxes of the 
PAH features and the [Nell] line, there is a clear progression in relative strength of the higher 
ionization lines [SIV] 10.5 /mi and [Nelll] 15.6 /im , from strongest in the GALEX sources, 
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intermediate in the Markarian sources, to weakest in the Spitzer sources. The [SIII] 18.7 
/im line is also strong in the Markarian sources compared to the Spitzer sources, although 
S/N is too low for measurement in the GALEX average spectrum. The average spectra 
indicate, therefore, that the ultraviolet-selected samples contain hotter ionizing stars than 
does the extreme infrared-selected Spitzer sample. 

The average spectra also show increased dust continuum for the GALEX sample com- 
pared to the other samples, easily seen by the higher continuum for wavelengths > 15 /im but 
also measurable by an increased continuum at ~ 10 /im. That dust emission increases while 
the hardness of the ionization also increases is consistent with having warmer dust within 
the HII region in the sources with harder ionization. Comparisons are plotted for individual 
sources in Figures 7 and 8. 

In Figure 7, the ratio of total flux in [Nelll] to PAH 11.3 /im is compared with EW of 
PAH 1 1.3 Atm (in /jm ) . Thi s plot includes all three samples plus the nearby starburst galaxies 
in the 



Brandl et al. 



Bernard-Salas et al. 



J2006h sample as measured with the high resolution spectroscopy in 



( 20091 ). [NeIII]/PAH is a measure of hardness of the ionizing radiation 



and EW(11.3 /im) shows the strength of the dust continuum relative to the strength of PAH. 
The median values of log ([NeIII]/PAH), including limits, are -0.29, -1.01, -1.13, and -1.25 
for the GALEX, Markarian, Spitzer, and nearby starbursts. Both Figure 5 and these results 
show decreasing EW(11.3 /im) PAH with increasing [NeIII]/PAH. The average EW(11.3 /im) 
for the Spitzer sample is 0.65 /mi, for the Markarian sample is 0.55 /im, and for the GALEX 
sample is 0.50 /an. 

The values for EW(11.3 /im) indicate that the Markarian starbursts have 18% more 
dust luminosity in the continuum at ~ 11 /im relative to PAH luminosity than has the 
Spitzer sample, and the GALEX sample has 30% more dust luminosity than the Spitzer 
sample. A similar result is found by measuring dust continuum directly. For example, 
the value of continuum f„(10 /im) in the normalized average spectra is 0.16/0.29/0.40 for 
Spztzer/Markarian/GALEX. The ratio of dust continuuum at 15 /im to PAH peak at 7.7 
/im is 0.48/0.68/0.91 for Spztzer/Markarian/GALEX. These results are all consistent with 
the simple interpretation that the harder ionizing radiation in the ultraviolet-selected star- 
bursts results in stronger dust emission at mid-infrared wavelengths, compared to the extreme 
Spitzer sample. 

Figure 8 compares the [NeIII]/PAH ratio with luminosity of the PAH features, measured 
by vL u (7.7 /im) in erg s _1 . There is no trend with luminosity for this ratio, but the result 
shows as in Figure 7 the systematically stronger [Nelll] compared to PAH in the UV-selected 
GALEX and Markarian samples compared to the Spitzer infrared-selected sample. 
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Obscuration in Starbursts 



Beca use the starbursts at highest redshift have b een discovered in the rest frame ul- 
traviolet ( IBouwens et al.ll2009l ; iReddy and Steidell 120091 ). understanding their intrinsic lumi- 
nosities requires corrections for dust obscuration. The mea sure of o b scura tion is the greatest 
uncertainty in measuring SFRs. ICalzetti et al.l ( 120071 ) and ICalzettil (120081 ) have summarized 
the various methods for determining such corrections and their application to measurement 
of SFRs. 

When observing the rest frame ultraviolet, the obscuration is estimated b y comparing 
the o bserved spectral slope (or "reddening") with an assumed intrinsic slope (e.g. lBurgarella et al 
20051 ). This is the source of the correction applied to the LBGs discovered at the highest 
redshifts. This has led to the result that the dust correction remains small for lower luminos- 
ity galaxies, which yields the con clusion that low lumin osity galaxies progressively dominate 
the luminosity density for z > 3 (IBouwens et al.ll2009n . 



Such obscuration corrections apply only for the portion of the ultraviolet luminosity 
that is able to emerge from the obscuration. Completely hidden sources cannot be measured 
in this way. This is the primary reason why infrared measures of starbursts are so important, 
because infrared luminosities are much less subject to extinction. In our discussion below, 
we present a new, direct comparison between the parameters most directly measuring the 
starburst luminosity as observed in the infrared (PAH luminosity) and in the ultraviolet (con- 
tinuum luminosity). These parameters are compared to determine empirically the relation 
for unobscured starbursts, and then to use this relation to determine obscuration corrections 
for obscured starbursts. 

The PAH emission and the ultraviolet monochromatic continuum provide two spectral 
"features" which are intrinsic to the same starburst. The f v (7.7 /iin) arises from the pho- 
todissociation region immediately surrounding the ongoing starburst and directly excited by 
photons from the starburst. This PAH emission is intrinsic to the photodissociation region 
and is present regardless of whether any dust obscures the starburst. It is just as intrin- 
sic to the starburst as radiation from the hot stars or emission lines from the HII region. 
Both f„(7.7 /im) and f^(153 nm) features are affected by the same obscuring dust in the 
surrounding cold molecular cloud, but the ultraviolet suffers heavy extinction whereas the 
PAH feature suffers negligible extinction. 

If we can determine empirically the intrinsic ratio between the PAH emission and ul- 
traviolet continuum in unobscured sources, we can then simply use the observed ratio in 
obscured sources to determine how much extinction applies to the ultraviolet. This empiri- 
cal determination is made using Figures 11, 12, and 14, discussed below, which illustrate the 
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most extreme, limiting values of f u (7.7 /^m)/f I/ (153 nm) found in any of the starbursts. In all 
three figures, this limit is log[f„(7.7 //m)/f iy (153 nm)] = 1. We will adopt this as the intrinsic 
value for sources with negligible obscuration of ^(153 nm). If these extreme sources also 
have some ultraviolet obscuration, then the final values of obscuration which we derive are 
lower limits, and sources are even more obscured in the ultraviolet than our results indicate. 



5.1. Comparison of Star Formation Rates 



The ultimate objective of studying the luminosities of starbursts is to determine the 
SFR in starbursts and to follow the resulting cosmic evolution of the SFR density. Various 
assumptions enter the transformation from the luminosity of a starburst into a quantita- 
tive measure of SFR in M Q yr _1 . For infrared measures based on luminosities of the dust 
continuum, the SFR is determined by equating the bolometric luminosity of the starburst 
to the luminosity absorbed and reradiated by dust. For ultraviolet luminosities, the SFR 
derives from the total luminosity of the hot stars in the starburst measured at a rest frame 
ultraviolet wavelength. 

For sources with a mid-infrared spectrum arising only from a starburst with no AGN 
contamination, the luminosity uL u (7.7 /im) relates empirically to the total i nfrared lumi 
nosity of the starburst, L ir , according to log L ir = log [z/L„(7.7 //m)] + 0.78 (IHouck et al. 
20071 ). This transformation has no dependence on luminosity, as is seen by combinin g the 
IRAS-derived L{ r for the local, low luminosity starburst sam ple of iBrandl et all ( 120061 ) with 
the high lum inosity, high reds hift submillimeter galaxies in iPope et all ( 120081 ) . Using the 
relation from Kennicutt ( 19981 ) to convert L ir into a SFR , this result for L ir yields that log 
[SFR(PAH)] = log [i/L„(7.7 /im)] - 42.5 7±0.2, for uLJ7.7 mm) i n ergs s" 1 and SFR in M 
yr _1 . This conversion was confirmed in ISargsyan and Weedmanl (120091 ) who found no offset 
between this measure of SFR and that derived from radio continuum observations of the 



same starbursts (jCondonlll992l ) 



SFR(UV) is deter mined using the rel ation SFR(UV)= 1.08xlO- 28 L„(153 nm), for 1^(153 
nm) in ergs s _1 Hz -1 (jSalim et all 120071 ). For our comparisons, SFR(UV) is the observed 
value without applying any extinction corrections to the ultraviolet. For comparison to 
SFR(PAH), we convert this to log [SFR(UV)] = log i/L„(FUV) - 43.26, for i/L„(FUV) 
in ergs s _1 and SFR in M yr _1 , using f^(153 nm) from Tables 1, 3, and 5. The ratio 
SFR(PAH)/SFR(UV) is a measure of the fraction of luminosity absorbed by dust and seen 
in the infrared compared to the fraction which escapes absorption and can be seen in the 
ultraviolet. 
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Figure 9 shows SFRs determined from 7.7 /mi PAH luminosity compared to SFRs deter- 
mined from observed GALEX UV continuum for the three samples of starbursts. The line is 



reprod uced form the full sample of 287 infrared-selected starbursts in lSargsyan and Weedman 



((20091). The GALEX and Markarian UV-selected samples show similar values of SFR(PAH) /SFR(UV), 
with median log SFR(PAH)/SFR(UV) = 0.8. If the discrepancy between SFR(PAH) and 
SFR(UV) is caused only by the obscuration of the UV, then this result indicates that only 
~ 15% of the intrinsic UV luminosity is observed, even for these ultraviolet-selected sam- 
ples. Much greater obscuration is found for the Spitzer-selected sample, with median log 
SFR(PAH)/SFR(UV) = 2.4, indicating that only 0.4% of the intrinsic ultraviolet luminos- 
ity emerges. The GALEX sample represents the most extreme UV-selected starbursts, and 
these are systematically more luminous than the Markarian sources. 

In Figure 10, we transform this comparison to use M(UV) as a measure of luminosity. 
M(UV) is an AB magnitude, deter mined with the fundamental definition m.AB = -2.51og i v 



48.60, for fj, in erg cm" ' 2 s" 1 Hz" 1 dQkelll974f ). M(UV) is determined as the m(UV) which 



arises using L^(153 nm) to determine f^(153 nm) at a distance of 10 pc. The line shown is 
the linear least squares fit to all of the points. It is notable that this line has the opposite 
slope of the line in Figure 9. This result indicates that objects selected using ultraviolet 
observations of luminosity select in favor of sources having less obscuration. As a result, 
the dependence of obscuration on luminosity has the opposite sign of the dependence found 
when starbursts are selected in the infrared. Figure 10 shows the similarity of GALEX and 
Markarian sources in terms of obscuration as measured by SFR(PAH)/SFR(UV), but shows 
that the GALEX sample is systematically more luminous. 



5.2. Comparing Infrared and Ultraviolet Starburst Luminosities to Determine 

Obscuration 

While the ultimate use of infrared or ultraviolet spectra of starbursts is to determine 
the star formation rates, the dependence of this result on the amount of obscuration means 
that it is useful to discuss only those parameters which can measure the obscuration. In 
this way, we hope to determine a quantitative measure of obscuration in starbursts which is 
independent of the many assumptions required to relate observed luminosity to SFR from 
either ultraviolet or infrared spectra. 

Therefore, we introduce a new comparison which is only between the observed param- 
eters used to determine PAH or UV luminosities. This avoids any of the uncertainties in 
converting to SFRs. These parameters are the observed f^(7.7 /mi) and fj,(153 nm) (Tables 
1, 3, and 5). We assume neglible extinction of the PAH, so that the PAH luminosity rep- 



- 16 - 



resents the in trinsic luminosity of the starburst. This is justified by the extinction curve of 



Draind ( 119891 ). which shows small PAH extinction compared to optical or UV extinction. For 
example, from Draine, A(Ha)/E(J-K) = 6.5, but A(7.7/im)/E(J-K) = 0.07, so 5 magnitudes 
of Ha extinction (1% emerges) corresponds to only 0.05 mag of 7.7 /im extinction (95% 
emerges) . 

Figure 11 shows the comparison of infrared PAH flux density to ultraviolet flux density, 
f„(7.7 / um)/f I/ (153 nm), as a function of luminosity uL u (7.7 /im). The distribution of points 
among the three samples illustrates an empirical result that the smallest value of log[f„(7.7 
/im)/f v (153 nm)] = 1. We take this as the empirical value which relates PAH luminosity to 
ultraviolet luminosity for a starburst with negligible obscuration of the UV. This is the most 
important result of the comparisons shown in Figure 11. It is an empirical result which can 
then be used to compare ultraviolet and infrared spectra for a quantitative determination of 
obscuration in the ultraviolet spectra. 

Figure 12 shows the same comparison but using M(UV) as the luminosity measure. This 
Figure also shows the empirical result that the limiting value of log[f !y (7.7 /im)/f v (153 nm)] 
= 1, so we take this to be the value for an unobscured source. The UV-selected sources 
in both Markarian and GALEX samples have similar values of log[f„(7.7 //m)/f !y (153 nm)], 
with a median of log[f^(7.7 //m)/f t/ (153 nm)] = 1.8, compared to the median for the Spitzer 
sample of 3.4. 

Figures 11 and 12 show the strong offset in the samples resulting from the selection 
criteria. The infrared-selected starbursts have median values of i P (7.7 /im)/f ;/ (153 nm) larger 
by about a factor of 100 compared to the ultraviolet-selected starbursts. This is a direct 
indication that the infrared-selected sample shows extreme obscuration in the ultraviolet. 
Demonstrating obscuration in this way avoids the uncertainties and assumptions required 
when comparing SFRs determined from the different samples. 

Using the limit for f u (7.7 /im)/f !y (153 nm) as representing sources without obscuration, 
such that ^[^(7.7 /im)/f v (153 nm)] = 1 means no obscuration, Figures 11 and 12 then 
give the correction factor which can be used to determine the obscuration correction for 
a source. The fraction of UV luminosity intrinsically emitted compared to the luminos- 
ity observed, UV(intrinsic)/UV(observed), is given by log[UV(intrinsic)/UV(observed)] = 
log[f„(7.7 /im)/f,(153 nm)] - 1. 
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5.3. Obscuration Corrections to M(UV) 

The correction for dust obscuration determined above in section 5.2 is compared with 
luminosity M(UV) in Figure 13. Figure 13 illustrates our most important results regarding 
the selection effects arising from obscuration by dust for discovering starbursts . The lumi- 
nosities of infrared-selected and ultraviolet-selected samples overlap, but the most luminous 
ultraviolet sources (the GALEX sample, diamonds in Figure 13) appear luminous because 
they have low obscuration. There is a very large range between the minimally obscured 
ultraviolet sources and the most obscured infrared sources, with the latter having dust cor- 
rection factors exceeding 100 at M(UV) = -18, compared to a correction of less than a factor 
of 10 at this M(UV) for the ultraviolet-selected sources. 

We show in Figure 13 the linear fit for the most obscured Spitzer starbursts (upper line), 
the sample in Table 5, but we also plot in Figure 13 the full sampl e of 287 infrared-selected 



starbu rsts (center line) having IRS and GALEX measures from ISargsyan and Weedman 



(120091 ). The distribution of these points shows that the ultraviolet-selected starbursts from 
the Markarian and GALEX samples (fit by the lower line) generally lie at the extremes of 
the full infrared-selected sample. 

The fits to the sources in Figure 13 give measures of obscuration which apply for star- 
bursts so obscured that most of the ultraviolet luminosity is hidden. These fi ts are indepen- 



dent o f mea sures derived only from reddening of ultraviolet observations, as in lBouwens et al. 



(120091 ) and iBurgarella et all (120051 ). Even in the ultraviolet-selected samples, the obscura- 
tion is significant, being a factor of 8 to 10 with little dependence on luminosity. The lin- 
ear least squares fit to the ultraviolet selected samples is logfUVf intrinsic) / UV (observed)] = 



0.07(±0.04)M(UV)+2.09±0.69. This obscuration is about twice that used bv lReddv and Steidel 



(120091 ) to correct for LBGs with z > 3, although the M(UV) in Figure 13 overlap with the 



M(UV) of these high redshift sources. 

For the full infrared-selected Spitzer sample (central line in Figure 13), The linear least 
squares fit is log[UV(intrinsic)/UV(observed)] = 0.17(±0.02)M(UV)+4.55(±0.4); the fit to 
the extreme infrared-selected Spitzer sample in Table 5 (upper line) is log [UV(intrinsic)/UV(observed)] 
= 0.14(±0.03)M(UV)+4.68±0.50. 

We make another estimate of effects of dust obscuration in sections 6 and 7, by deter- 
mining the relation between M(UV) and total bolometric luminosities L ir . This comparison 
is possible because we have an empirical scaling between Lj r and uL u (7.7 /mi). 
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5.4. Ratio Ha/H/3 as Measure of Obscuration compared to 

Infrared/Ultraviolet 

A method often used to determine obscuration corrections for starbursts is the measure 



of reddening from the Ha/H/3 ratio (ICalzetti et all 120071 ). This method works only for 



the emission lines which are able to emerge from the starbursts and so cannot account for 
totally obscured starbursts. It is useful, therefore, to attempt a comparison between the 
obscuration which would be deduced from Ha/H/3 and that deduced from f v (7.7 //m)/f I/ (153 
nm). Our primary objective of this comparison is to obtain another measure of the value 
fj/(7.7 /im)/f^(153 nm) for an unobscured starburst. 

Figure 14 compares Ha/H/3 determined from SDSS spectra of the starbursts (Tables 1, 
3, 5) with the ratio f„(7.7 /^m)/f„(153 nm) discussed in section 5.2. The relative Ha/H/3 line 
fluxes were measured using AfA taking fA at the peak of the line in the SDSS spectra. Values 
given in Tables 1, 3, and 5 are for all sources with SDSS spectra available. The intrinsic 
Ha/H/3 ratio in an unobscured source should be 2.87, which indeed agrees with lowest value 
in Figure 14. 

The comparisons in Figure 14 confirm the determination of log[f^(7.7 / um)/f l/ (153 nm)] 
= 1 for an unobscured source, found above in section 5.2. The line is a linear least squares 
fit to all of the points. The result that there is a trend line shows that there is indeed 
some relation between reddening observed for the Balmer lines and the obscuration shown 
by comparing PAH and UV luminosities. 

To quantify this relation, we show in Figure 15 a measure of obscuration for Ha com- 
pared with the obscuration deduced from comparing f u (7.7 /^m)/f„(153 nm) discussed in sec- 



tion 5.2. From ICalzetti et all (120071 ) . k(Ha) = 2.47 and K(H/3)-K(Ha)=1.16. Knowing the 
intrinsic value Ha/H/3 = 2.87, this differential extinction means that 0.471og Ha(intrinsic) = 
1.471og Ha(observed)-log H/3(observed)-0.46. From this, we derive a value Ha(intrinsic) /Ha (observed). 

This fraction of escaping Ha is compared to the fraction of escaping UV by using the 
empirical result in section 5.2 that log[UV(intrinsic)/UV(observed)] = log[f„(7.7 //m)/f I/ (153 
nm)]-l. The comparison is shown in Figure 15. 

Figure 15 shows the expected result for highly obscured regions that the obscuration 
of the intrinsic UV is much greater than the obscuration which would be deduced from the 
Balmer lines. This is because the Balmer lines which can be observed emerge from outside 
of the most obscured regions, whereas our measure of ultraviolet obscuration applies to the 
total obscuration of the most hidden regions. 
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5.5. Comparing Starburst and Stellar Luminosities for Starburst Galaxies 



Understanding the formation of galaxies requires, in addition to a measure of the 
SFR, a determination of the mass of stars already present. Such masses can be estimated 
from observations at various rest-fra me wavelengths if s t ellar models and mixtures ar e as- 
sumed, as explained, fo r example, in iReddy and Steidell (120091 ). lLonsdale et al.l (120091 ) and 
Bussmann et al.l (120091 ) . For the starburst samples in the present paper, we compare the 
luminosities of starbursts with the luminosities of the underlying stellar component for the 
galaxy. 

This comparison is done by showing in Figure 16 the ratio of f v (7.7 fim) at the peak 
of the PAH feature to the rest frame f„(1.6 fJ>m), which arises from the stellar continuum. 
For the Markarian, GALEX, and Spitzer samples, the stellar luminosity is measured from 
the Two Micron All Sky Survey (2MASS; Skrutskie et al. 2006) H band flux densities; no 
corrections for redshift are made because there are only small differences between the H 
flux density and K (2.2 /^m) flux density (Tables 1, 3, and 5.) Figure 16 is also useful for 
comparison of PAH luminosities among the various samples, including the DOGs at z ~ 2. 

In Figure 16, we include numerous starburst DOGs from the optically faint, high redshift 
samples selected at f„(24 /im). All have published IRS spectra that allow the measure of 
f„(7.7 jum), and have photometry with the Spitzer Infrared Array Camera (IRAC; Fazio et 
al. 2004) that allows mea s urement of rest fr a me 1.6 /im luminosi ty. T hese sources are fo und 



in 



Weedman et al.l (I2006h : lYan et all (120071 ): iFarrah et al.l (120081 ) and iDesai et all (120091 ) 



These sources represent the most luminous sources known as meas ured with i/L„(7.7 /im) 
so th ey show the most luminous, most obscured, starburst ULIRGs (jWeedman and Houck 
20081 ). Most of these sources were selected as "bump" sources, for which an excess continuum 
is observed at 4. 6 um to 5.8 /im from giant and supergiant stars having lower opacity at rest 
frame ~ 1.6 /im (jSimpson and Eisenhardtlll999l ). The comparison in Figure 16 is with IRAC 
5.8 /im fluxes to estimate rest frame f„(1.6 /im) because fluxes at 4.6 /im and 5.8 /im are so 
similar that no redshift corrections are made. 

Because the high redshift starbursts were selected based both on f^(24 /im) and ^(5.8 
/zm), it is expected that selection effects will restrict the range of f v (7.7 //m)/f^(1.6 jum), the 
measure of starburst luminosity /stellar luminosity. No such selection effect can apply for the 
lower redshift samples from the present paper, however, because 2MASS magnitudes played 
no role in choosing these sources. 

The similarity in the starburst to stellar ratio is striking among all of the samples. The 
range of 10 4 in PAH luminosity compares to a range of only 4 relative to stellar luminosity. 
We do not attempt to calculate stellar masses because of uncertainties in converting a 1.6 
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/j,m stellar luminosity into a stellar mass. More importantly, the homogeneity of the f u (7.7 
/im)/f^(1.6 /im) ratio among the samples in Figure 16 raises a fundamental question of what 
such a stellar mass measure at 1.6 /im means. We wonder if the stellar luminosity at 1.6 
jttm arises from supergiant stars very recently evolved from the starburst itself and does not 
reflect the underlying stellar mass of old stars. 

There is some evidence for this interpretation if we compare the ratios i v {7.7 /^m)/f ;y (1.6 
jum) among the different samples. This evidence is that the median value of the ratio 
f„(7.7 jim)/i y {l.Q jum) scales with obscuration determined from i v {7.7 /im)/f !y (153 nm), as 
measured in section 5.2. The ultraviolet-selected GALEX and Markarian sample combined 
have a median i v (7.7 /im)/f„(1.6 /im) = 8 (counting limits), the low redshift infrared-selected 
extreme Spitzer sample in Table 5 (local DOGs) has median f u (7.7 /im)/fj,(1.6 jum) = 13, 
and the high redshift DOGs sample of the most luminous, most dusty starbursts have median 
f„(7.7 /im)/t(1.6 H = 20. 

While extinction at 1.6 /mi is much less than in the ultraviolet, some extinction can also 
be expected at 1.6 /xm for th e most obscured sources. Because extinction at 7.7 /im is less 



than at 1.6 fim (jDrainelll989l ). the ratio i u (7.7 /im)/f J/ (1.6 /im) should scale with the amount 
of obscuration, more obscured sources having a higher ratio, if the 1.6 /im stellar luminosity 
arises within the starburst and is also obscured by dust. This is the trend which is observed 
among the samples, from most obscured Spitzer sources to least obscured ultraviolet sources. 
By contrast, we would not expect such a scaling if the 1.6 /xm luminosity arises from an 
underlying stellar component which is outside of the obscured starbursts. 



6. Comparing Bolometric Luminosities L ir with M(UV) 



The fundamental question for measurement of SF Rs is which parameter best represent s 



the integrated SFR over the life of the starburst (e.g. iCalzettil 120081 ; iRosa-GonzalezI |2002| ) . 
This is the estimate required to determine the rate at which galaxies form. Different param- 
eters are sensitive to different stars in the starbursts and to starbursts of different ages. For 
example, the Ha emission depends on the luminosity of the most massive, short lived stars 
with the greatest ionizing radiation shortward of the Lyman limit, but L^(153 nm) measures 
luminosity from hot stars longward of the Lyman limit, so can be influenced by lower mass, 
longer lived stars. Comparing these measures would not give the same result for starbursts 
with different mass distributions for the ionizing stars or different ages for the starburst. 



The preferred parameter is the total luminosity of all stars in the starburst, which does 
not depend on what fraction of the stellar luminosity ionizes hydrogen. This parameter is 
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measured by Lj r , the total bolometric luminosity over all wavelengths, and it is dominated 
by infrared luminosity because most starburst luminosity is ab sorbed by obscu ring dust. 
This is the basis of measures of SFR using infrared luminosities ( IKennicuttl 119981 ). 



The measure of SFR that w e have used from uL v (7.7 /xm) is related empirically to Lj r 
by using the starburst galaxies in lBrandl et al.l (120061 ). which all have IRAS flux measures al- 
lowin g the determinatio n of L ir . The transformation adopted is log L ir = log[z/L Jy (7.7/im)] + 
0.78 (IHouck et al.ll2007l ). and t his is confirmed t o extend to starbursts at z > 2 using submil- 
limeter determinations of L ir (IPope et al.ll2008l ). (The L ir in Brandl e t al. is a total 8-1000 
/ im in frared luminosity as estimated from IRAS fluxes according to ISanders and Mirabel 
(119961 ).) Using this transformation, we determine log Lj r compared to M(UV) for all sources 
in our samples. 

Results are shown in Figures 17 and 18, comparing the ultraviolet-selected samples to 
the infrared-selected samples. In Figure 17, the lines which are fit compare the extreme 
Spitzer sample of Table 5 to the ultraviolet Markarian+GALEX samples, the extremes 
of infrared selection and ultraviolet selection for starbursts. This Figure demonstrates in 
alternative fashion to the discussion in section 5.2 that there exists a population of very 
obscured starbursts which have larger luminosities Lj r than found among ultraviolet-selected 
starbursts and which have faint M(UV) because of this obscuration. 

The offset between the fits in Figure 17 indicates that ultraviolet-selected starbursts of 
a given M(UV) can underestimate bolometric luminosities by about a factor of 30 over all 
luminosities when compared with the most obscured starbursts. Another way to state the 
result is that starbursts of the same L ir can be found at all -16 > M(UV) > -21, and those 
which appear faintest in M(UV) are so because of obscuration. If there are large populations 
of starbursts at z > 3 which are as dusty and obscured as these extreme Spitzer sources, 
the conclusion that only low luminosity galaxies exist at such redshifts is incorrect. 

In Figure 18, we show the sample s compared with the fit found for the full infrared- 
selected sample of 287 starbursts from ISargsyan and Weedmanl (120091 ) . In this case, differ- 
ences between lines for the ultraviolet-selected sample (lower line) and the infrared-selected 
sample (upper line) are not as large as in Figure 17. The offset gives another estimate of 
obscuration corrections for ultraviolet-selected galaxies when compared to infrared-selected 
samples for comparison with the obscuration corrections in section 5.2. 

For the samples and fits in Figures 17 and 18, we conclude the following. For ultraviolet- 
selected galaxies, log L ir = -(0.33±0.04)M(UV)+4.52±0.69. For the full Spitzer infrared- 
selected starburst sample, log L ir = -(0.23±0.02)M(UV)+6.99±0.41. For the extreme infrared- 
selected Spitzer sample in Tables 5 and 6, log L ir = -(0.25±0.03)M(UV)+7.2±0.48, for L ir 
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in L and M(UV) the AB magnitude at rest frame 153 nm. 

All of these interpretations depend on an accurate determination of L ir . A crucial 
requirement is to improve the calibration of L ir to ^L J/ (7.7/im). We can make one test of 
the adopted calibration using our new IRS spectra of Markarian galaxies, because most have 
IRAS fluxes. We show in Figure 19 the result of this ratio for both the Brandl starbursts 
and for our new observations of Markarian galaxies in Table 2. 

Results in Figure 19 show that, within the uncertainties, the ultraviolet-selected Markar- 
ian galaxies have the same transformation from Lj r to zAL , y (7.7//m) as we have explained and 
adop ted in previous analyses for Spitzer-selected sources ( IHouck et al. IWeedman and Houckl 



20081 ). This is an important confirmation that this transformation and the subsequent use 
of L ir to determine SFRs is not modified by starbursts with the hotter ionizing continuum 
seen in the Markarian starbursts. Nevertheless, improvements for this calibration using far 
infrared and millimeter observations of starbursts with a variety of luminosities and redshifts 
are essential. 



7. Obscuration Corrections and Evolution for Starburst Galaxies at High 

Redshift 



Our motivation for comparing ultraviolet-selected and infrared-selected starbursts is to 
determine an independent measure of obscuration by dust, because the value of such a dust 
correction is crucial for measuring the SFR density. Present conclusions (summarized in 
Bouwens et al. 2009), which track the SFR density to z ~ 7, indicate that obscuration for 
luminous galaxies, M(UV) < -20, is substantial, requiring a correction factor to ultraviolet 
luminosities of 6±2.5. The overall dust correction factor to total SFR density is taken as 
only a factor of ~ 2, however, because the correction is smaller for low luminosity galaxies 
which dominate the observed SFR density. As concluded by Bouwens et al., "Because these 
same lower luminosity galaxies dominate the luminosity density in the UV continuum, the 
overall dust extinction correction remains modest at all redshifts and the evolution of this 
correction with redshift is only modest." 

Having a correct measurement of dust obscuration is crucial to determine the true SFR 
density as a function of redshift because the form of evolution for SFR density depends on 
this correction. For exam ple, the observed decrease in the SFR density from z ~ 2 to z ~ 7 in 
Reddy and Steidell (120091 ) is comparable to the factors applied for dust correction for the most 
luminous sources. Many assumptions and e xplanations ente r in th e det ermination of such 
dust corrections, as thoroughly discussed in iBouwens et al.l ( 120091 ) and iReddy and Steidel 
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( 120091 ). Fundamentally, the adopted corrections depend on interpreting as reddening by 
dust the observed differences in the slopes of the ultraviolet continuum for LBGs of different 
luminosities. 

The ultimate measure of SFR derives from the corrected values of Lj r which are finally 
adopted after applying the corrections. The consequences of different methods for deter- 
mining corrections can be compared, therefore, simply by comparing the derived L ir as a 
function of ultraviolet luminosity, M(UV). These are the comparisons we now make using 
our methodology from section 6 for determining L ir based on infrared PAH luminosities. 

Results are shown in Figure 20. We consider that Figure 20 summarizes the most 
important results of our present paper. The relation betwee n L ir and M(TJV) for t he various 



samples which we consider is compared to that adopted in iBouwens et al.l (120091 ) for LBGs 
at high redshift used to determine the evolution of SFR density for z > 2. In Figure 20, the 
line with asterisks is taken from Figure 10 of Bouwens et al. for z = 2.5. The comparison 
of this line with the points in Figure 20 indicates that the result adopted for high redshift 
LBGs really applies only to the lower envelope for all starbursts. 

Differences between this line and the other lines in Figure 20 show differences that would 
also arise for resulting measures of SFR density. The adopted results for LBGs agree with 
the full infrared-selected sample at M(UV) ~ -21, but differ by a factor > 10 at M(UV) = 
-17. Differences are greater at fainter M(UV) for reasons discussed in section 5.2, that the 
faintest sources in the ultraviolet are those with the greatest obscuration by dust. 

The crucial question that results from Figure 20 is: which line represents the true 
population of starburst galaxies at high redshift? Are the dust properties determined for our 
local samples (z < 0.5) similar to properties in high redshift starbursts of similar M(UV)? 
Answering this question for z > 2 is not feasible at present, because there are no significant 
samples of sources discovered at such redshifts using their dust emission. 

A direct comparison could be made at z ~ 2 because of the large population of DOGs 
discovered by Spitzer at about this redshift. A strong selection effect leads to a median 
redshift z = 1.8±0.2 for the starburst DOGs with f I ,(24 /im) > 0.5 mJy, using the high 
redshift, dusty starbursts shown in Figure 16 (asterisks). This selection effect arises when 
the strong PAH 7.7 fim feature falls within the MIPS 24 /mi band used for the surveys 
that find these DOGs. When the 24 jum selection is also coupled to a criterion of having 
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This selection using "bump" sources means that existing Spitzer surveys covering tens of 
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deg 2 , such as the First Lo ok Survey (|Fadda et al. 1120061) . the Bootes survey (IDey et al.ll2008l ) 



and the SWIRE surveys (ILonsdale et al. 112004] ). can yield reliable lists of dusty starbursts in 



a well defined redshift interval near z ~ 2 , even if spectroscopic redshifts are not available for 
all of the sources. There is also a straightforward empirical transformation in this redshift 
interval that f^(7.7 /an) = (1.5±0. 3)^(24 /mi), for DOGs like those in Figure 16. The space 
density and luminosity function of such starbursts could be determined to f^(24 /mi) > 0.3 
mJy, for which the f u (7.7 /mi) ~ 0.5 mJy, or log vL u (7.7 /mi) = 45.2 at z = 1.8. 

From Figures 17 and 18, this value of vL v (7 .7 /mi) corresponds to -16 > M(UV) > 
-22 , depending on whether the source is highly obscured or has lower obscuration such as 
in the ultraviolet-selected samples. This range of M(UV) is easily accessible in surveys for 
UVLGs at z ~ 2, so such surveys in the same regions where are available the Spitzer surveys 
would make possible a direct comparison between infrared and ultraviolet selection of the 
starburst population at z ~ 2. This would provide a definitive answer regarding obscuration 
corrections and measures of true SFR densities as a function of M(UV). 

Even though existing wide-area infrared surveys reach only the top end of the bolometric 
luminosity function at high redshifts, it is encouraging that tracking the most luminous 
sources gives results for the form of evolution over all redshifts that are consistent with results 
found from total SFR densities which reach deeper into luminosity functi ons. We reach this 
conclu sion by comparing the evolution of the SFR density as illustrated in lReddy and Steidel 
(120091 ). their Figure 10, with the luminosity evolution factor of (T+z) 2,5 found from Spitzer 
surveys for the most luminous starbursts (jWeedman and Houck 2008 ). This scaling for 
luminosity evolution mat ches, within the uncertain ties, the scaling of total SFR density for 
< z < 2 summarized in Reddy and Steidel ( 2009 ). 



8. Summary and Conclusions 

1. We have selected starburst galaxies discovered using both infrared and ultraviolet 
selection for comparison of spectral parameters measured with the Spitzer Infrared Spectro- 
graph to ultraviolet photometry from GALEX. Our primary objective was to determine how 
the luminosity and dust obscuration of starbursts depend on selection technique. Sources 
have z < 0.5 and cover a luminosity range of ~ 10 4 ; the highest luminosities are comparable 
to starburst luminosities found at z ~ 3. 

2. We measure and illustrate spectra of 26 Markarian galaxies, 23 ultraviolet luminous 
galaxies discovered with GALEX, and 50 galaxies with the most extreme infrared/ultraviolet 
ratios discovered with Spitzer. We find (section 4.3) that these samples differ in the nature 
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of the ionizing stars, with hotter stars in the ultraviolet-selected sources as evidenced by 
stronger [Nelll] emission and stronger mid-infrared dust continuum (section 4, Figures 6 - 
8.) 

3. It is found that a strong selection effect arises for the ultraviolet-selected samples: the 
most luminous UV sources appear luminous because they have the least obscuration, not be- 
cause they have the largest bolometric luminosities (section 5.3). Comparisons between rest- 
frame infrared fluxes f u (7.7 //m) arising from PAH emission and ultraviolet fluxes fj,(153 nm) 
arising from the stellar continuum allow an empirical method for determining obscuration 
in starbursts and dependence of this obscuration on infrared or ultraviolet luminosity. Even 
sources selected in the ultraviolet and which resemble Lyman Break Galaxies require obscura- 
tion corrections of about a factor of 10. These results are contrary to those adopted in studies 
of LBGs, which find that the most luminous galaxies in M(UV) are also the most obscured. 
Obscuration correction for the ultraviolet-selected Markarian+GALEX sample has the form 
log[UV(intrinsic)/UV(observed)] = 0.07(±0.04)M(UV)+2.09±0.69; the full infrared-selected 
Spitzer sample is log[UV(intrinsic)/UV(observed)] = 0.17(±0.02)M(UV)+4.55(±0.4); and is 
log [UV(intrinsic)/UV(observed)] = 0.14(±0.03)M(UV)+4.68±0.50 for the extreme infrared- 
selected Spitzer sample (Figure 13). 

4. Analysis of obscuration determined from the Hct/H/3 ratio in SDSS spectra confirms 
the trends observed by comparing infrared/ultraviolet ratios, although obscuration measured 
by Hct/H/3 is much less than measured by infrared/ultraviolet (section 5.4 and Figure 15). 

5. Stellar luminosities measured at rest frame 1.6 ^m are exceptionally uniform com- 
pared to starburst luminosity over a redshift range < z < 2 and luminosity range of 10 4 , but 
the ratio from different samples indicates that the 1.6 /iin stellar luminosity is also obscured 
by dust in proportion to the ultraviolet obscuration. This leads to our suggesting that stellar 
luminosities measured at 1.6 /im reveal the recently evolved stars from the starburst rather 
than an underlying old population (section 5.5 and Figure 16). 

6. We examine the relation of total bolometric luminosity Li r to M(UV) for infrared- 
selected and ultraviolet-selected samples. Spectroscopy of Markarian galaxies confirms the 
transformation assumed for vL u (7.7 /im) to L ir and implies that the hotter ionizing stars 
in the uv-selected samples do not affect this calibration. For ultraviolet-selected galaxies, 
log L ir = -(0.33±0.04)M(UV)+4.52±0.69. For the full infrared-selected sample, log L ir = - 
(0.23±0.02)M(UV)+6.99±0.41, and log L ir = -(0.25±0.03)M(UV)+7.2±0.48 for the extreme 
infrared-selected Spitzer sample, all for L ir in L Q and M(UV) the AB magnitude at rest frame 
153 nm (section 6 and Figure 18). 

7. Our most important summary result (section 7 and Figure 20) indicates that the 
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bolometric luminosities deduced for starburst galaxies using LBGs at z > 2 are insufficient, 
by a factor of at least fO for M(UV) ~ -17. This also implies that obscuration corrections 
by factors of two to three applied for 3 < z < 7 to track the formation and evolution of 
the earliest galaxies are insufficient, if dusty galaxies exist at such epochs comparable to 
those with z < 3. We propose an observational test which can determine more accurate dust 
corrections and L ir for large samples of infrared-selected starbursts and LBGs at z < 2. This 
test is possible because existing Spitzer surveys allow reliable discovery of dusty starbursts 
with z = 1.8 ± 0.2 based only on photometric properties. 
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Table 1. Spectroscopic Data for Markarian Sample 



No. 


Name 


J2000 coordinates 


z a 


EW(6.2) b 


U7.7) c 


f(11.3) d 


EW(11.3) 


f(15.5) c 


f„(24) f 


f„(153)* 


Ha/H/3) h 


f(Hy 


f(K) 










fim 


mJy 




/an 




mJy 


^Jy 




mJy 


mJy 


1 


Mkn 353 


010316.58+222032.6 


0.014 


0.71 


176(2.04) 


23.18 


0.78 


1.47 


422 


266 




18.1 


17.0 


2 


Mkn 1158 


013459.36+350222.1 


0.015 


0.69 


47.3(1.19) 


3.43 


0.50 


1.27 


160 


403 




5.6 


5.3 


3 


Mkn 363 


015058.56+215948.0 


0.010 


0.65 


64.5(1.41) 


6.64 


0.70 


1.16 


145 


772 




9.4 


8.0 


4 


Mkn 1012 


015657.73-052411.7 


0.020 


0.48 


90.8(1.41) 


9.35 


0.54 


<0.52 


466 


218 




14.6 


14.7 


5 


Mkn 603 


030856.85-025718.4 


0.007 


0.61 


258(1.63) 


26.67 


0.48 


13.88 


1790 


1140 




17.5 


18.9 


6 


Mkn 14 


081059.31+724740.3 


0.008 


0.59 


28.1(1.47) 


2.62 


0.51 


0.51 


120 


884 




5.9 


4.7 


7 


Mkn 87 


082140.17+735917.5 


0.009 


0.52 


48.7(1.28) 


4.68 


0.74 


0.07 


71.3 


187 




14.6 


12.3 


8 


Mkn 90 


083000.10+524148.9 


0.014 


0.52 


30.3(1.32) 


2.88 


0.65 


0.16 


80.4 


194 


5.5 


4.8 


5.0 


9 


Mkn 93 k 


083642.35+661358.3 


0.017 


0.36 


48.7 ••• 


2.77 


0.19 


3.02 


367 


147 




3.5 


3.4 


10 


Mkn 105 


092026.34+712416.0 


0.012 


0.62 


13.9 ••• 


0.88 


0.50 


0.17 


45.3 


120 




2.9 


2.2 


11 


Mkn 21 


094830.73+575814.4 


0.027 


0.45 


24.4(1.54) 


2.66 


0.70 


0.1 


40.3 


305 


7.2 


7.3 


o&.l 


12 


Mkn 140 


101628.23+451919.4 


0.004 


0.51 


13.9 ••• 


1.02 


0.58 


0.60 


43.1 


738 


3.6 


2.4 


o 

,1.9 


13 


Mkn 31 


101942.83+572524.9 


0.026 


0.50 


22.8 • • • 


1.74 


0.66 


0.15 


58.0 


135 


4.6 


4.3 


3.8 


14 


Mkn 149 k 


103801.78+641559.0 


0.005 


0.24 


17.0 ••• 


1.53 


0.23 


<0.08 


165 


163 


5.3 


7.3 


7.4 


15 


Mkn 165 


111835.60+631641.0 


0.011 


0.61 


16.8(1.57) 


1.95 


0.78 


0.26 


23.0 


785 


4.9 


8.1 


6.6 


16 


Mkn 193 k 


115528.50+573944.0 


0.021 


<0.04 


<2.9 • • • 


0.07 


0.03 


0.4 


60.4 


332 




<0.3 




17 


TV /T1 'Tp y|k 

Mkn 764 


121600.04+124114.0 


0.066 


0.47 


32.9(1.49) 


3.23 


0.64 


0.07 


42.9 


242 


7.2 


3.6 


3.6 


18 


Mkn 439 


122436.30+392259.0 


0.005 


0.61 


132.0(1.80) 


17.86 


0.72 


0.71 


409.0 


3760 




22.6 


18.8 


19 


Mkn 769 


122525.54+162811.8 


0.006 


0.59 


194(1.10) 


14.48 


0.61 


5.62 


638 


5770 




29.1 


21.7 


20 


Mkn 799 


140045.70+591942.0 


0.009 


0.51 


207(1.67) 


22.54 


0.51 


0.46 


857 




6.9 


31.6 


43.6 


21 


Mkn 286 


141926.97+713517.4 


0.026 


0.64 


187(1.46) 


16.16 


0.62 


1.14 


597 


1220 




14.5 


13.0 


22 


Mkn 480 


150633.73+423829.2 


0.019 


0.66 


101(1.45) 


8.91 


0.72 


0.29 


170 


523 


5.2 


8.5 


9.2 


23 


Mkn 489 


154421.68+410509.8 


0.032 


0.62 


77.0(1.57) 


7.7 


0.58 


1.14 


217 


972 


4.3 


3.6 


3.2 


24 


Mkn 691 


154658.90+175308.0 


0.011 


0.65 


106(1.34) 


8.13 


0.52 


0.88 


355 




4.4 


12.7 


11.2 


25 


Mkn 863 k 


155625.90+090318.8 


0.043 


0.30 


47.6(1.16) 


3.02 


0.38 


0.14 


103 


341 


5.3 


4.8 


6.5 


26 


NGC 7714 1 


233614.10+020919.0 


0.010 


0.59 


354(1.17) 


25.35 


0.36 


7.36 


2150 


4190 




21.8 


20.4 



a Redshift from PAH features in IRS spectra. 

b Rcst frame equivalent width of 6.2 ^im PAH feature fit with single gaussian on a linear continuum within rest wavelength range 5.5 /im to 6.9 /zm . 

c Flux density in mJy at peak of 7.7 (im PAH feature. The value in parenthesis is the "stitching" factor by which the SL1 spectrum has been multiplied 
to match the LL2 spectrum to account for a slightly resolved source that gives larger LL2 flux because the SL1 slit is 3.7 " wide, and the LL2 slit is 
10.5 " wide. 

d Total flux of PAH 11.3 /jm in units of 1CT 20 W enr 2 s" 1 . 
e Total flux of [Nelll] 15.5 /um in units of 1(T 20 W cnT 2 s" 1 . 
f Continuum flux density at rest frame 24 /im. 

s Flux density in /jjy observed at wavelength (l+z)153 nm interpolated with a power law between flux densities in the GALEX catalog measured 
with the FUV filter (observed 134-179 nm) and the NUV filter (observed 177-283 nm). 

h Observed emission line ratio Ha/H/3 measured by us from SDSS spectra using ratios of AfA for i\ at the line peak shown in the SDSS spectra. 

'Observed flux density in the H band of 2MASS using conversion that zero magnitude corresponds to 1024 Jy. 

1 Observed flux density in the K band of 2MASS using conversion that zero magnitude corresponds to 667 Jy. 

k SDSS spectrum shows classification as AGN, as discussed in section 4.2. 

'All spectra from our new Spitzer program 50834, except NGC 7714 from archival program 14. 



Table 2. Luminosities and Star Formation Rates for Markarian Sample 



No. 


Name 


J2000 coordinates 


logL(7.7/zm) a 


\oeUUV) h 


M(UV) C 


log SFR(PAH) d 


log SFR(UV) e 


L,r f 








lot? prff s — ^ 


lot? ertr s — ^ 




loffM^ vr 


lOEfMr^ VT 


lot? erff s 


1 


Mkn 353 


010316.58+222032.6 


43.46 


42.34 


-16.03 


0.89 


-0.92 


44.07,(0.70) 


2 


Mkn 1158 


013459.36+350222.1 


42.95 


42.58 


-16.63 


0.38 


-0.68 


<43.65,(0.70) 


3 


Mkn 363 


015058.56+215948.0 


42.73 


42.51 


-16.45 


0.16 


-0.75 


<43.24,(1.20) 


4 


Mkn 1012 


015657.73-052411.7 


43.48 


42.57 


-16.61 


0.91 


-0.69 


44.41,(0.83) 


5 


Mkn 603 


030856.85-025718.4 


43.02 


42.37 


-16.09 


0.45 


-0.89 


43.99,(0.79) 


6 


Mkn 14 


081059.31+724740.3 


42.17 


42.37 


-16.11 


-0.40 


-0.89 


<43.04,(0.77) 


7 


Mkn 87 


082140.17+735917.5 


42.52 


41.80 


-14.68 


-0.05 


-1.46 


<43.23,(0.61) 


8 


Mkn 90 


083000.10+524148.9 


42.70 


42.20 


-15.69 


0.13 


-1.06 


<43.41,(0.49) 


9 


Mkn 93 s 


083642.35+661358.3 


43.07 


42.25 


-15.80 






43.85,(0.75) 


10 


Mkn 105 


092026.34+712416.0 


42.22 


41.86 


-14.83 


-0.35 


-1.40 


<43.01,(0.53) 


11 


Mkn 21 


094830.73+575814.4 


43.17 


42.97 


-17.61 


0.60 


-0.29 


<43.67,(0.36) 


12 


Mkn 140 


101628.23+451919.4 


41.26 


41.69 


-14.40 


-1.31 


-1.57 


<41.86,(0.39) 


13 


Mkn 31 


101942.83+572524.9 


43.11 


42.59 


-16.64 


0.54 


-0.67 




14 


Mkn 149 s 


103801.78+641559.0 


41.55 


41.23 


-13.25 






<42.46,(0.82) 


15 


Mkn 165 


111835.60+631641.0 


42.23 


42.60 


-16.68 


-0.34 


-0.66 


<42.57,(0.29) 


16 


Mkn 193 s 


115528.50+573944.0 


<42.03 


42.79 


-16.73 






43.50,(0.57) 


17 


Mkn 764 s 


121600.04+124114.0 


44.09 


43.66 


-19.32 






<44.30,(0.23) 


18 


Mkn 439 


122436.30+392259.0 


42.44 


42.60 


-16.66 


-0.13 


-0.66 


43.28,(0.58) 


19 


Mkn 769 


122525.54+162811.8 


42.76 


42.90 


-17.43 


0.19 


-0.36 


43.54,(0.60) 


20 


Mkn 799 


140045.70+591942.0 


43.15 






0.58 




43.99,(0.50) 


21 


Mkn 286 


141926.97+713517.4 


44.03 


43.54 


-19.04 


1.46 


0.28 


44.71,(0.77) 


22 


Mkn 480 


150633.73+423829.2 


43.49 


42.90 


-17.43 


0.92 


-0.36 


44.07,(0.83) 


23 


Mkn 489 


154421.68+410509.8 


43.82 


43.63 


-19.24 


1.25 


0.37 


<44.65,(0.91) 


24 


Mkn 691 


154658.90+175308.0 


43.03 






0.46 




43.85,(0.64) 


25 


Mkn 863 s 


155625.90+090318.8 


43.87 


43.43 


-18.75 






44.26,(0.04) 


26 


NGC 7714 


233614.10+020919.0 


43.44 


43.22 


-18.22 


0.87 


-0.04 


44.19,(0.71) 



a Rest frame luminosity of 7.7 /im PAH feature in erg s _1 determined as L(7.7 fim) = v\j v (7 .7 /im) = 47rDi 2 [^/(H-z)]/ J/ (7.7 
/im), for v corresp onding t o 7.7 /xm and using fj,(7.7 /an) at observed wavelength of the 7.7 yum feature from Table 1. Luminosity 
distances D l from IWrightl(l2006h : |http:// www. astro . ucla.edu / ~ wright / CosmoCalc . htrnl] for H = 71 kms 1 Mpc 1 , Om=0.27 and 
Q A =0.73. (Log [i/L„(L )] = log [i/L„(erg s" 1 )] - 33.59.) 



b Rest frame ultraviolet continuum luminosity at 153 nm in erg s 1 determined as L(UV) 
Table at observed wavelength (l+z)153 nm. 



^^(153 nm), using ^(153 nm) from 



C M(UV) is an absolute AB magnitude at co ntinuum w avelength 153 nm, determined with the fundamental definition vciab = 
-2.51ogf I/ - 48.60, for i v in erg cm~ 2 s _1 Hz -1 (|Okelll974l ). M(UV) is determined as the m(UV) which arises using L^(153 /*m) to 
determine f„(153 nm) at a distance of 10 pc. 

d Star formation rates in units of M Q yr _1 determined as log [SFR(PAH)] = log L(7.7 /im) - 42.57. 

c Star formation rates in units of M© yr -1 determined as log [SFR(UV)] = log L(UV) - 43.26. 

f Total infrared luminosity determined from IRAS fl ux densities by taking t otal in frared flux from 8-1000 ^m to be 1.8[13.48S(12/im) 
+ 5.16S(25^m) + 2.58S(60^m) + S(100^m)] from Zanders and Mirabell |l996h . where S is IRAS flux density in Jy. Value in 
parenthesis is the ratio f„(IRS 25 /zm)/f„ (IRAS 25 /im) determined by using synthetic photometry with SMART for the IRAS 25 
/im band to determine (IRS 25 /im) from the IRS spectrum in Figure 1. 



S SDSS spectrum shows classification as AGN, as discussed in section 4.2. 



Table 3. Spectroscopic Data for GALEX Sample 



No. 


NcLlTLG 


J2000 coordinates 




z a 


EW(6.2) b 


f,X7 7Y 


ffll 3) d 


EW(11.3) 


f(15.5) c 


f„(24) f 


i,(153) g 


Ha/H/3) h 


i(H)' 1 


i(K\ 












/zm 


m Jv 




fun 




m Jv 


11,1 V 




m Jv 


m Jv 


1 


SDSS k 


001054.84+001451.3 


0. 


.244 


0.15 


6.3 


0.34 


0.42 


<0.01 


13.5 


26.5 


8.21 


0.8 


1.1 


2 


SDSS 


005439.79+155446.9 





.239 


0.27 


3.3 


0.15 


0.33 


<0.01 


7.1 


13.2 


5.68 


<0.3 




3 


SDSS k 


005527.46-002148.8 


0. 


.168 


0.08 


6.6 


0.20 


0.10 


0.25 


50.3 


77.9 


3.94 


0.3 


0.4 


4 


SDSS k 


015028.41+130858.4 





.147 


0.16 


18.1 


0.34 


0.26 


0.31 


93.4 


82.0 


4.01 


0.5 


0.5 


5 


SDSS 


032846.00+011150.8 





.143 


0.46 


1.2 


0.09 


0.65 


<0.03 


6.1 


31.3 


3.76 


<0.3 




6 


SDSS 


035733.99-053719.7 





.204 


0.42 


1.9 


0.07 


0.37 


<0.06 


7.2 


12.6 


4.14 


<0.3 




7 


SDSS 


040208.87-050642.1 


0. 


.138 


0.74 


0.5 


0.03 


0.35 


<0.02 


1.4 


45.2 


3.88 


<0.3 




8 


SDSS 


080232.34+391552.6 


0. 


.267 


0.66 


4.1 


0.16 


0.43 


0.01 


12.3 


34.8 


5.84 


<0.3 




9 


SDSS 


080844.27+394852.3 


0. 


.093 


0.53 


9.2 


0.52 


0.29 


<0.02 


51 


160.5 


4.69 


0.9 


0.6 


10 


SDSS k 


082001.72+505039.1 





.219 


0.04 


3.6 


0.03 


0.06 


0.05 


23.9 


30.5 


3.70 


<0.3 




11 


SDSS 


082550.95+411710.2 





.156 


0.46 


1.7 


0.09 


0.50 


0.05 


4.7 


42.2 


4.06 


0.2 


0.4 


12 


SDSS k 


083803.72+445900.3 





.145 


<0.09 


0.4 


0.03 


0.40 


0.03 


3.6 


<77.0 


3.09 


<0.3 




13 


SDSS 


092159.38+450912.3 





.236 


0.23 


11.3 


0.42 


0.51 


0.06 


44.7 


78.3 


5.37 


0.4 


0.6 


14 


SDSS 


092336.45+544839.2 


0. 


.227 


0.58 


1.1 


0.08 


0.64 


0.03 


4.1 


44.5 


4.19 


<0.3 




15 


SDSS k 


092600.41+442736.1 





.182 


<0.17 


<0.3 


<0.002 


<0.01 


0.06 


5.9 


<94.1 


3.24 


<0.3 




16 


SDSS 


093813.50+542825.1 





.101 


0.18 


2.3 


0.11 


0.16 


0.29 


28.3 


249.8 


3.38 


0.3 


0.5 


17 


SDSS 


102613.97+484458.9 


0. 


.157 


0.91 


1.3 


0.09 


0.57 


0.07 


5.9 


85.8 


3.34 


<0.3 




18 


SDSS 


124819.75+662142.7 


0. 


.260 


0.67 


1.4 


0.08 


0.52 


0.04 


5.5 


33.8 


4.08 


<0.3 




19 


SDSS 


135355.90+664800.6 


0. 


.200 


0.90 


3.0 


0.20 


0.81 


0.05 


8.2 


100.0 


3.80 


0.4 


0.6 


20 


SDSS 


143417.15+020742.5 





.181 


0.61 


3.7 


0.24 


0.68 


0.06 


9 


50.7 


4.47 


<0.3 




21 


SDSS 


210358.75-072802.5 


0. 


.136 


0.38 


21.1 


0.11 


0.19 


0.17 


128.9 


84.8 


5.20 


1.7 


1.9 


22 


SDSS 


214500.25+011157.5 





.205 


0.66 


3.9 


0.22 


0.77 


0.08 


8.5 


32.7 


4.13 


<0.3 




23 


SDSS k 


231812.99-004126.1 





,253 


<0.03 


<6.4 


0.17 


0.31 


0.17 


31.4 


<67.6 


3.68 


<0.3 





a Rcdshift from PAH features in IRS spectra. 

b Rest frame equivalent width of 6.2 /zm PAH feature fit with single gaussian on a linear continuum within rest wavelength range 5.5 /zm to 
6.9 /zm . 

c Flux density in mJy at peak of 7.7 ^m PAH feature. 
d Total flux of PAH 11.3 ^m in units of KT 20 W cnT 2 s" 1 . 
c Total flux of [Nelll] 15.5 fim in units of 1CT 20 W cm" 2 s" 1 . 
'Continuum flux density at rest frame 24 /zm. 

s Flux density in /zJy observed at wavelength (l+z)153 nm interpolated with a power law between flux densities in the GALEX catalog 
measured with the FUV filter (observed 134-179 nm) and the NUV filter (observed 177-283 nm). 

h Observed emission line ratio Ha/H/3 measured by us from SDSS spectra using ratios of Af>, for f\ at the line peak shown in the SDSS 
spectra. 

'Observed flux density in the H band of 2MASS using conversion that zero magnitude corresponds to 1024 Jy. 
J Observed flux density in the K band of 2MASS using conversion that zero magnitude corresponds to 667 Jy. 
k SDSS spectrum or IRS spectrum shows classification as AGN, as discussed in section 4.2. 



Table 4. Luminosities and Star Formation Rates for GALEX Sample 



No. 


Name 


J2000 coordinates 


logL(7.7/zm) a 


logLfUV) b 


M(UV) C 


log SFR(PAH) d 


log SFR(UV) 








lot? orec s 1 


off frer s ~^ 




loe']VI+ vr 


loffM^ vr ^ 


1 


SDSS f 


001054.84+001451.3 


44.54 


43.87 


-19.84 






2 


SDSS 


005439.79+155446.9 


44.24 


43.55 


-19.04 






3 


SDSS f 


005527.46-002148.8 


44.22 


44.00 


-20.17 


1.65 


0.74 


4 


SDSS f 


015028.41+130858.4 


44.54 


43.90 


-19.93 






5 


SDSS 


032846.00+011150.8 


43.35 


43.46 


-18.82 


0.78 


0.20 


6 


SDSS 


035733.99-053719.7 


43.86 


43.38 


-18.64 


1.29 


0.12 


7 


SDSS 


040208.87-050642.1 


42.92 


43.59 


-19.14 


0.35 


0.33 


8 


SDSS 


080232.34+391552.6 


44.43 


44.07 


-20.35 


1.86 


0.81 


9 


SDSS 


080844.27+394852.3 


43.84 


43.78 


-19.64 


1.27 


0.52 


10 


SDSS f 


082001.72+505039.1 


44.20 


43.83 


-19.75 






11 


SDSS 


082550.95+411710.2 


43.57 


43.67 


-19.34 


1.00 


0.41 


12 


SDSS f 


083803.72+445900.3 


42.87 


43.86 


-19.83 






13 


SDSS 


092159.38+450912.3 


44.76 


44.31 


-20.94 


2.19 


1.05 


14 


SDSS 


092336.45+544839.2 


43.72 


44.03 


-20.24 


1.15 


0.77 


15 


SDSS f 


092600.41+442736.1 


<43.0 


44.15 


-20.56 






16 


SDSS 


093813.50+542825.1 


43.31 


44.05 


-20.30 


0.74 


0.79 


17 


SDSS 


102613.97+484458.9 


43.47 


43.98 


-20.13 


0.90 


0.72 


18 


SDSS 


124819.75+662142.7 


43.93 


44.03 


-20.25 


1.36 


0.77 


19 


SDSS 


135355.90+664800.6 


44.04 


44.26 


-20.84 


1.47 


1.00 


20 


SDSS 


143417.15+020742.5 


44.04 


43.88 


-19.87 


1.47 


0.62 


21 


SDSS 


210358.75-072802.5 


44.54 


43.85 


-19.79 


1.97 


0.59 


22 


SDSS 


214500.25+011157.5 


44.18 


43.80 


-19.68 


1.61 


0.54 


23 


SDSS f 


231812.99-004126.1 


<44.58 


44.31 


-20.94 







a Rest frame luminosity of 7.7 /im PAH feature in erg s" 1 determined as L(7.7 /jm) = i/L v (7.7 
/im) — 47rD^ 2 [^/(l+z)]/ l/ (7.7 /im), for ^ corresponding to 7.7 /im and using f„(7.7 /i m) at o bserve d 
wavelength of the 7.7 /im feature from Table 1. Luminosity distances D L from Wright ( 20061) : 



http://www.astro.ucla.edu/~wright/CosmoCalc.htmf for H = 71 kms 1 Mpc l , f^/=0.27 and £7a=0.73. 



(Log [uL v (L @ )] = log [i/L„(erg s" 1 )] - 33.59.) 

b Rcst frame ultraviolet continuum luminosity at 153 nm in erg s _1 determined as L(153 nm) = zTj„(153 
nm), using f„(153 nm) from Table 3 at observed wavelength (l+z)153 nm. 

C M(UV) is an absolute AB magnitude at continuum wavelength 153 nm, de termined with the fundamental 
definition tclab = -2.51ogf y - 48.60, for f„ in erg cm~ 2 s _1 Hz -1 |Okelll974r) . M(UV) is determined as the 
m(UV) which arises using L„(153 /im) to determine f„(153 nm) at a distance of 10 pc. 

d Star formation rates in units of M Q yr _1 determined as log [SFR(PAH)] = log L(7.7 /im) - 42.57. 

c Star formation rates in units of M Q yr _1 determined as log [SFR(UV)] = log L(UV) - 43.26. 

f SDSS spectrum or IRS spectrum shows classification as AGN, as discussed in section 4.2. 



Table 5. Spectroscopic Data for Spitzer Sample 



No. Name J2000 coordinates z a EW(6.2) b i„(7.7) c f(11.3) d EW(11.3) f(15.5)° f„(24) f f„(153) s Ha/H/3) h f(i?)' i(Ky 













/iiri 


mJy 




fim 




mJy 






mJy 


mJy 


1 


SST24 


021849.78-052158.4 





292 


0.71 


6.2 


0.24 


0.53 


<0.01 


11.6 


3.5 




<0.3 




2 


SST24 


022205.05-050536.7 





259 


0.52 


7.7 


0.35 


0.81 


<0.02 


12.7 


1.9 




0.5 


0.6 


3 


SST24 


022345.09-054234.4 





142 


0.55 


10.7 


0.6 


0.76 


<0.04 


13.7 


8.2 




0.8 


0.6 


4 


SST24 


022447.00-040851.4 





096 


0.36 


12.5 


0.6 


0.74 


<0.02 


8.7 


3.4 




1.8 


1.9 


5 


SST24 


022548.25-050051.8 





145 


0.42 


4.2 


0.2 


0.55 


<0.01 


14.3 


0.5 




0.4 


0.5 


6 


SST24 


022549.79-040025.1 





043 


0.36 


23.8 


1.41 


0.50 


0.27 


69.6 


15.1 




2.4 


2.3 


7 


SST24 


022600.02-050145.1 





203 


0.55 


6.4 


0.32 


0.72 


<0.06 


9.6 


3.0 




0.6 


0.6 


8 


SST24 


022637.83-035841.8 





073 


0.38 


5.1 


0.28 


0.52 


<0.02 


18.5 


3.2 




0.7 


0.7 


9 


SST24 


022738.55-044702.9 





171 


0.51 


5.8 


0.33 


0.68 


<0.03 


11.9 


3.8 




0.4 


0.5 


10 


SST24 


104016.34+570845.9 





115 


0.64 


7.3 


0.35 


0.70 


<0.02 


7.2 


2.3 




0.5 


0.6 


11 


SST24 


104729.91+572843.3 





230 


0.64 


11.2 


0.59 


0.79 


<0.02 


10.5 


6.0 


7.71 


0.6 


0.8 


12 


SST24 


104907.19+565715.1 





072 


0.75 


16.7 


0.83 


0.76 


<0.03 


11.7 


14.1 


6.19 


1.0 


1.0 


13 


SST24 


105336.93+580350.7 





455 


0.39 


4 


0.13 


0.22 


<0.03 


6.2 


3.0 




<0.3 




14 


SST24 


110133.82+575205.8 





282 


0.54 


8.2 


0.33 


0.58 


<0.01 


14.2 


1.3 




0.5 


0.5 


15 


SST24 


142504.04+345013.7 





078 


0.59 


16.6 


0.9 


0.78 


0.04 


14.5 


10.8 


6.44 


1.6 


1.4 


16 


SST24 


142554.57+344603.2 





035 


0.56 


38.6 


2.19 


0.72 


0.11 


60.2 


27.5 


6.03 


3.5 


3.3 


17 


SST24 


143039.27+352351.0 





089 


0.47 


14.4 


0.71 


0.53 


<0.06 


49.1 


6.1 


10.58 


1.2 


1.4 


18 


SST24 


143119.79+353418.1 





033 


0.68 


39 


2.39 


0.78 


0.18 


39.7 


42.3 


5.82 


3.4 


3.3 


19 


SST24 


143126.81+344517.9 





084 


0.54 


13.7 


0.76 


0.72 


<0.01 


14.1 


3.3 


7.58 


1.2 


1.1 


20 


SST24 


160038.84+551018.4 





143 


0.42 


9.7 


0.49 


0.43 


<0.04 


30.6 


2.2 




0.6 


0.8 


21 


SST24 


160341.36+552612.5 





145 


0.59 


8.9 


0.52 


0.75 


0.03 


8.0 


4.3 




0.9 


1.2 


22 


SST24 


160358.21+555504.9 





333 


0.39 


3.9 


0.04 


0.08 


0.03 


11.5 


0.3 




<0.3 




23 


SST24 


160408.22+542531.2 





264 


0.54 


7 


0.3 


0.75 


0.03 


9.4 


0.6 




<0.3 




24 


SST24 


160408.36+545812.7 





064 


0.52 


20.1 


1.13 


0.67 


<0.07 


32.8 


5.4 


9.99 


2.6 


2.6 


25 


SST24 


160832.64+552927.0 





065 


0.59 


8.7 


0.43 


0.69 


<0.01 


7.4 


1.9 




0.7 


0.9 


2G 


SST24 


160907.60+552428.3 





064 


0.53 


9.6 


0.46 


0.76 


<0.06 


9.6 


8.6 




0.7 


0.8 



Table 5 — Continued 



No. 


Name 


J2000 coordinates 




z a 


EW(6.2) b 
/iin 


f,(7.7) c 
mJy 


f(11.3) d 


EW(11.3) 
fim 


f(15.5)° 


f,(24) f 
mJy 


1.(153)6 
,ujy 


Ka/Rf3) h 


W 
mJy 


i(K) 
mJy 


27 


SST24 


160908.34+552241.2 





086 


0.57 


11.6 


0.72 


0.92 


0.11 


8.4 


6.4 




1.0 


1.1 


28 


SST24 


160937.59+541259.6 





085 


0.65 


12 


0.69 


0.81 


<0.02 


6.4 


7.1 




1.0 


1.5 


29 


SST24 


161823.14+552721.3 





083 


0.64 


36.3 


1.91 


0.77 


0.08 


36.2 


31.4 




1.7 


1.8 


30 


SST24 


161827.74+552209.2 





082 


0.53 


8.5 


0.39 


0.73 


<0.05 


14.6 


0.8 




0.7 


0.8 


31 


SST24 


161950.56+543715.4 





144 


0.54 


10.3 


0.51 


0.68 


<0.03 


10.2 


5.2 




0.7 


0.9 


32 


SST24 


162038.14+553521.6 





190 


0.57 


9.4 


0.47 


0.86 


<0.02 


14.7 


6.8 




0.5 


0.6 


33 


SST24 


162210.98+550253.5 





030 


0.42 


31.4 


1.42 


0.62 


0.08 


55.7 


3.7 




4.7 


4.4 


34 


SST24 


163001.55+410953.3 





123 


0.58 


7.5 


0.39 


0.73 


<0.04 


11.0 


1.1 




<0.3 




35 


SST24 


163359.13+405305.0 





031 


0.58 


15 


0.82 


0.79 


0.07 


12.0 


13.7 


5.86 


1.0 


1.2 


36 


SST24 


163506.08+411038.2 





078 


0.52 


14.8 


0.59 


0.63 


<0.03 


18.7 


1.1 


9.84 


0.9 


1.0 


37 


SST24 


163608.15+410507.8 





167 


0.41 


14.6 


0.51 


0.59 


0.1 


26.5 


10.5 


9.53 


0.8 


1.0 


38 


SST24 


163715.60+414934.2 





122 


0.52 


5.8 


0.33 


0.64 


<0.05 


12.2 


4.6 


6.99 


0.8 


0.9 


39 


SST24 


163751.95+401503.6 





072 


0.69 


12.4 


0.63 


0.83 


<0.03 


10.8 


5.9 




0.7 


0.8 


40 


SST24 


163808.55+403214.2 





225 


0.46 


9.2 


0.37 


0.57 


0.04 


26.0 


3.4 




0.3 


0.5 


41 


SST24 


164211.93+410816.4 





143 


0.51 


9 


0.43 


0.62 


0.06 


20.9 


4.8 




0.5 


0.6 


42 


SST24 


171316.55+583235.5 





080 


0.52 


9.5 


0.5 


0.74 


0.02 


9.4 


9.7 


9.20 


1.6 


1.6 


43 


SST24 


171414.89+585221.9 





166 


0.59 


7.1 


0.32 


0.65 


0.04 


14.9 


1.4 




0.4 


0.6 


44 


SST24 


171446.54+593400.3 





130 


0.62 


14.7 


0.7 


0.79 


0.05 


10.4 


3.3 




0.8 


1.1 


45 


SST24 


171542.00+591657.4 





119 


0.54 


12 


0.5 


0.57 


0.1 


43.9 


6.6 




1.0 


1.3 


46 


SST24 


171614.52+595423.4 





153 


0.58 


10.4 


0.57 


0.72 


<0.02 


13.2 


9.4 


6.76 


0.7 


0.6 


47 


SST24 


171711.17+602709.4 





110 


0.53 


6.8 


0.36 


0.66 


<0.04 


12.6 


1.3 




0.6 


0.7 


48 


SST24 


171933.40+592743.2 





142 


0.45 


9.1 


0.37 


0.61 


0.03 


13.3 


2.4 




0.6 


0.7 


49 


SST24 


172043.33+584026.5 





125 


0.54 


9.7 


0.57 


0.73 


<0.04 


15.0 


5.7 


7.51 


1.1 


1.1 


50 


SST24 


172400.61+590228.2 





179 


0.44 


6.2 


0.29 


0.45 


<0.03 


23.7 


6.7 




0.5 


0.6 



a Redshift from PAH features in IRS spectra. 

b Rest frame equivalent width of 6.2 ^m PAH feature fit with single gaussian on a linear continuum within rest wavelength range 5.5 ^m to 
6.9 /im . 

c Flux density in mjy at peak of 7.7 /im PAH feature. 
d Total flux of PAH 11.3 ^m in units of 1(T 20 W cm" 2 s" 1 . 
°Total flux of [Nelll] 15.5 ^m in units of 1CT 20 W cm" 2 s" 1 . 
f Continuum flux density at rest frame 24 /im. 

s Flux density in /xJy observed at wavelength (l+z)153 nm interpolated with a power law between flux densities in the GALEX catalog 
measured with the FUV filter (observed 134-179 nm) and the NUV filter (observed 177-283 nm). 

h Observed emission line ratio Ha/H/3 measured by us from SDSS spectra using ratios of Af^ for i\ at the line peak shown in the SDSS 
spectra. 

'Observed flux density in the H band of 2MASS using conversion that zero magnitude corresponds to 1024 Jy. 
j Observed flux density in the K band of 2MASS using conversion that zero magnitude corresponds to 667 Jy. 



Table 6. Luminosities and Star Formation Rates for Spitzer Sample 



No. 


Name 


J2000 coordinates 


logL(7.7^m) a 


loeLCUV) 13 


M(UV) C 


log SFR(PAH) d 


log SFR(UV) 
loffM^ vr ^ 








lot? fT£? S 






lo£?M+ vr 


1 


SST24 


021849.78-052158.4 


44.70 


43.16 


-18.08 


2.13 


-0.10 


2 


SST24 


022205.05-050536.7 


44.68 


42.78 


-17.13 


2.11 


-0.48 


3 


SST24 


022345.09-054234.4 


44.29 


42.87 


-17.35 


1.72 


-0.39 


4 


SST24 


022447.00-040851.4 


44.00 


42.13 


-15.50 


1.43 


-1.13 


5 


SST24 


022548.25-050051.8 


43.90 


41.70 


-14.43 


1.33 


-1.56 


6 


SST24 


022549.79-040025.1 


43.57 


42.08 


-15.38 


1.00 


-1.18 


7 


SST24 


022600.02-050145.1 


44.38 


42.76 


-17.08 


1.81 


-0.50 


8 


SST24 


022637.83-035841.8 


43.37 


41.87 


-14.85 


0.80 


-1.39 


9 


SST24 


022738.55-044702.9 


44.19 


42.70 


-16.93 


1.62 


-0.56 


10 


SST24 


104016.34+570845.9 


43.93 


42.13 


-15.50 


1.36 


-1.13 


11 


SST24 


104729.91+572843.3 


44.74 


43.17 


-18.10 


2.17 


-0.09 


12 


SST24 


104907.19+565715.1 


43.87 


42.50 


-16.43 


1.30 


-0.76 


13 


SST24 


105336.93+580350.7 


44.91 


43.48 


-18.88 


2.34 


0.22 


14 


SST24 


110133.82+575205.8 


44.78 


42.69 


-16.90 


2.21 


-0.57 


15 


SST24 


142504.04+345013.7 


43.96 


42.46 


-16.33 


1.39 


-0.80 


16 


SST24 


142554.57+344603.2 


43.62 


42.16 


-15.58 


1.05 


-1.10 


17 


SST24 


143039.27+352351.0 


44.01 


42.32 


-15.98 


1.44 


-0.94 


18 


SST24 


143119.79+353418.1 


43.58 


42.30 


-15.93 


1.01 


-0.96 


19 


SST24 


143126.81+344517.9 


43.94 


42.00 


-15.18 


1.37 


-1.26 


20 


SST24 


160038.84+551018.4 


44.24 


42.31 


-15.95 


1.67 


-0.95 


21 


SST24 


160341.36+552612.5 


44.22 


42.61 


-16.70 


1.65 


-0.65 


22 


SST24 


160358.21+555504.9 


44.61 


42.27 


-15.85 


2.04 


-0.99 


23 


SST24 


160408.22+542531.2 


44.66 


42.29 


-15.90 


2.09 


-0.97 


24 


SST24 


160408.36+545812.7 


43.85 


41.98 


-15.13 


1.28 


-1.28 


25 


SST24 


160832.64+552927.0 


43.50 


41.54 


-14.03 


0.93 


-1.72 


26 


SST24 


160907.60+552428.3 


43.53 


41.82 


-14.73 


0.96 


-1.44 


27 


SST24 


160908.34+552241.2 


43.87 


42.31 


-15.95 


1.30 


-0.95 



Table 6 — Continued 



No. 


Name 


J2000 coordinates 


logL(7.7/mi) a 


logL(UV) b 


M(UV) C 


log SFR(PAH) d 


log SFR(UV) C 








log erg s _1 


log erg s -1 




logM Q yr _1 


logM Q yr _1 


28 


SST24 


160937.59+541259.6 


43.88 


42.35 


-16.05 


1.31 


-0.91 


29 


SST24 


161823.14+552721.3 


44.34 


42.98 


-17.63 


1.77 


-0.28 


30 


SST24 


161827.74+552209.2 


43.70 


41.37 


-13.60 


1.13 


-1.89 


31 


SST24 


161950.56+543715.4 


44.28 


42.69 


-16.90 


1.71 


-0.57 


32 


SST24 


162038.14+553521.6 


44.49 


43.05 


-17.80 


1.92 


-0.21 


33 


SST24 


162210.98+550253.5 


43.38 


41.15 


-13.05 


0.81 


-2.11 


34 


SST24 


163001.55+410953.3 


44.00 


41.87 


-14.85 


1.43 


-1.39 


35 


SST24 


163359.13+405305.0 


43.08 


41.75 


-14.55 


0.51 


-1.51 


36 


SST24 


163506.08+411038.2 


43.89 


41.46 


-13.83 


1.32 


-1.80 


37 


SST24 


163608.15+410507.8 


44.56 


43.12 


-17.98 


1.99 


-0.14 


38 


SST24 


163715.60+414934.2 


43.88 


42.49 


-16.40 


1.31 


-0.77 


39 


SST24 


163751.95+401503.6 


43.74 


42.12 


-15.48 


1.17 


-1.14 


40 


SST24 


163808.55+403214.2 


44.63 


42.91 


-17.45 


2.06 


-0.35 


41 


SST24 


164211.93+410816.4 


44.22 


42.65 


-16.80 


1.65 


-0.61 


42 


bbI24 


1 ti 01 ^? p"p" i roooor c 

171316.55+583235.5 


43.72 


42.43 


-16.25 


1.15 


-0.83 


43 


SST24 


171414.89+585221.9 


44.24 


42.24 


-15.78 


1.67 


-1.02 


44 


SST24 


171446.54+593400.3 


44.34 


42.40 


-16.18 


1.77 


-0.86 


45 


SST24 


171542.00+591657.4 


44.17 


42.62 


-16.73 


1.60 


-0.64 


46 


SST24 


171614.52+595423.4 


44.34 


43.00 


-17.68 


1.77 


-0.26 


47 


SST24 


171711.17+602709.4 


43.86 


41.86 


-14.83 


1.29 


-1.40 


48 


SST24 


171933.40+592743.2 


44.21 


42.34 


-16.03 


1.64 


-0.92 


49 


SST24 


172043.33+584026.5 


44.13 


42.60 


-16.68 


1.56 


-0.66 


50 


SST24 


172400.61+590228.2 


44.25 


42.99 


-17.65 


1.68 


-0.27 



a Rest frame luminosity of 7.7 /um PAH feature in erg s 1 determined as L(7.7 /im) = vh u (7.7 
/Ltm) = 47rD£ 2 [^/(l+z)]/iy(7.7 M m )i f° r v corresponding to 7.7 /im and using f„(7.7 /im) at observed 



wavelength of the 7.7 fan feature from Table 1. Luminosity distances from IWrightl toodl : 

http://www.astro.ucla.edu/~wright/CosmoCalc.html for H = 71 kms -1 Mpc -1 , fij^=0.27 and f2A=0.73. 
(Log [j/L„(L )] = log [vK(erg s" 1 )] - 33.59.) 

b Rest frame ultraviolet continuum luminosity at 153 nm in erg s _1 determined as L(153 nm) = i/L„(153 nm), 
using ^(153 nm) from Table 3 at observed wavelength (l+z)153 nm. 

C M(UV) is an absolute AB magnitude at continuum wavele ngth 153 n m, determined with the fundamental 
definition tclab = -2.51ogf !/ - 48.60, for i v in erg cm -2 s -1 Hz -1 (|Okelll974r ). M(UV) is determined as the m(UV) 
which arises using L^(153 fim) to determine f„(153 nm) at a distance of 10 pc. 

d Star formation rates in units of M Q yr _1 determined as log [SFR(PAH)] = log L(7.7 fim) - 42.57. 

°Star formation rates in units of M Q yr _1 determined as log [SFR(UV)] = log L(UV) - 43.26. 
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Fig. 1A. — Spectra of Markarian starbursts from Table 1. Markarian AGN are in Figure 4. 
All spectra are normalized to f u (7.7 /im) = 1.0 mJy, but zero points are displaced by 1 mJy 
for illustration. The zero flux level for each spectrum is located 1 mJy below the value at 
7.7 /xm. Spectra are labeled by running numbers in Table 1. Spectra are arranged in order 
of relative continuum strength at 24 /im, i.e. in order of ratio f„(24 /tm)/f I/ (7.7 /tm) in Table 
1. 
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wavelength (/xm) 

Fig. IB. — Spectra of Markarian starbursts from Table 1. Markarian AGN are in Figure 4. 
All spectra are normalized to i u (7.7 /im) = 1.0 mJy, but zero points are displaced by 1 mJy 
for illustration. The zero flux level for each spectrum is located 1 mJy below the value at 
7.7 /xm. Spectra are labeled by running numbers in Table 1. Spectra are arranged in order 
of relative continuum strength at 24 /im, i.e. in order of ratio i u (2A /tm)/f I/ (7.7 /xm) in Table 
1. 



-46- 



1 5 




Fig. 2A.— Spectra of GALEX starbursts from Table 3. GALEX AGN are in Figure 4. All 
spectra are normalized to f u (7.7 /im) = 1.0 mJy, but zero points are displaced by 1 mJy for 
illustration. The zero flux level for each spectrum is located 1 mJy below the value at 7.7 
fim. Spectra are labeled by running numbers in Table 3. Spectra are arranged in order of 
relative continuum strength at 24 /an, i.e. in order of ratio f^(24 //m)/f„(7.7 /im) in Table 3. 
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Fig. 2B.— Spectra of GALEX starbursts from Table 3. GALEX AGN are in Figure 4. All 
spectra are normalized to f^(7.7 /im) = 1.0 mJy, but zero points are displaced by 1 mJy for 
illustration. The zero flux level for each spectrum is located 1 mJy below the value at 7.7 
fim. Spectra are labeled by running numbers in Table 3. Spectra are arranged in order of 
relative continuum strength at 24 /an, i.e. in order of ratio f^(24 //m)/f„(7.7 /im) in Table 3. 
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Fig. 3 A. — Spectra of Spitzer starbursts from Table 5. All spectra are normalized to ^(7.7 
//m) = 1.0 mJy, but zero points are displaced by 1 mJy for illustration. The zero flux level for 
each spectrum is located 1 mJy below the value at 7.7 /an. Spectra are labeled by running 
numbers in Table 5. Spectra are arranged in order of relative continuum strength at 24 ^m, 
i.e. in order of ratio f^(24 /zm)/f„(7.7 ^m) in Table 5. 
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Fig. 3B. — Spectra of Spitzer starbursts from Table 5. All spectra are normalized to U(7.7 
/im) = 1.0 mJy, but zero points are displaced by 1 mJy for illustration. The zero flux level for 
each spectrum is located 1 mJy below the value at 7.7 /an. Spectra are labeled by running 
numbers in Tables 5. Spectra are arranged in order of relative continuum strength at 24 ^m, 
i.e. in order of ratio f^(24 /zm)/f„(7.7 ^m) in Table 5. 
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Fig. 3C. — Spectra of Spitzer starbursts from Table 5. All spectra are normalized to U(7.7 
/im) = 1.0 mJy, but zero points are displaced by 1 mJy for illustration. The zero flux level for 
each spectrum is located 1 mJy below the value at 7.7 /an. Spectra are labeled by running 
numbers in Table 5. Spectra are arranged in order of relative continuum strength at 24 ^m, 
i.e. in order of ratio f^(24 /zm)/f„(7.7 ^m) in Table 5. 
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Fig. 3D. — Spectra of Spitzer starbursts from Table 5. All spectra are normalized to ^(7.7 
//m) = 1.0 mJy, but zero points are displaced by 1 mJy for illustration. The zero flux level for 
each spectrum is located 1 mJy below the value at 7.7 /an. Spectra are labeled by running 
numbers in Table 5. Spectra are arranged in order of relative continuum strength at 24 ^m, 
i.e. in order of ratio f^(24 /zm)/f„(7.7 //m) in Table 5. 
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Fig. 3E. — Spectra of Spitzer starbursts from Table 5. All spectra are normalized to U(7.7 
/im) = 1.0 mJy, but zero points are displaced by 1 mJy for illustration. The zero flux level for 
each spectrum is located 1 mJy below the value at 7.7 /an. Spectra are labeled by running 
numbers in Table 5. Spectra are arranged in order of relative continuum strength at 24 ^m, 
i.e. in order of ratio f^(24 /zm)/f„(7.7 //m) in Table 5. 
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Fig. 4A. — Spectra of Markarian or GALEX AGN from Tables 1 and 3. All spectra are 
normalized to f u (7.7 /im) = 1.0 mJy, but zero points are displaced by 2 mJy for illustration. 
The zero flux level for each spectrum is located 1 mJy below the value at 7.7 /im. Spectra 
are labeled by running numbers in Tables 1 or 3. Spectra are arranged in order of relative 
continuum strength at 24 /im, i.e. in order of ratio f^(24 yum)/f l/ (7.7 /im) in Tables 1 or 3. 
Classification of these AGN is discussed in section 4.2. 
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Fig. 4B. — Spectra of Markarian or GALEX AGN from Tables 1 and 3. All spectra are 
normalized to f u (7.7 /im) = 1.0 mJy, but zero points are displaced by 2 mJy for illustration. 
The zero flux level for each spectrum is located 1 mJy below the value at 7.7 /im. Spectra 
are labeled by running numbers in Tables 1 or 3. Spectra are arranged in order of relative 
continuum strength at 24 /im, i.e. in order of ratio f^(24 yum)/f l/ (7.7 /im) in Tables 1 or 3. 
Classification of these AGN is discussed in section 4.2. 
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Fig. 5. — Top to bottom are average spectra for starbursts in the Markarian sample (Figure 
1), the GALEX sample (Figure 2), and the Spitzer sample (Figure 3). All spectra are 
normalized to f u (7.7 /im) = 1.0 mJy, but zero points are displaced by 1 mJy for illustration. 
The zero flux level for each spectrum is located 1 mJy below the value at 7.7 /im. For the 
Markarian sample, all sources are included in the average except the 5 AGN noted in Tables 
1 and 2. For the GALEX sample, all sources are included in the average except the 7 AGN 
noted in Tables 3 and 4. For the Spitzer sample, all sources are included because all are 
classified as starbursts from both SDSS and IRS spectra. 
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Fig. 6. — PAH 6.2 /im EW (in //m) compared with luminosity z/L^(7.7 /im) in erg s _1 
(log [z/Lv(L Q )] = log [z/L^erg s -1 )] - 33.59). Crosses are Spitzer starbursts from Figure 3. 
Triangles are Markarian starbursts from Figure 1, and circles are Markarian AGN from Figure 
4. Diamonds are GALEX starbursts from Figure 2, and squares are GALEX AGN from 
Figure 4. Small symbols show limits on undetectable PAH for Markarian and GALEX AGN. 
The Markarian starbursts are all found to satisfy the criterion of EW(6.2 /xm ) > 0.4 /j,m . 
The 3 GALEX sources with starburst classification but EW(6.2 /jm ) < 0.4 are especially 
interesting and are discussed in section 4.2. The Spitzer sample was initially chosen to select 
only sources with EW(6.2 /im ) > 0.4 ^m because of the well determined correlation between 
optical classification as starbursts and EW(6.2 //m ) > 0.4 jj,m . The Markarian starburst 
sample independently confirms this EW criterion, because the Balzano(1983) selection for 
the Markarian starbursts was defined using only optical criteria for starbursts. 



-57- 



=1 



< 
CL 



CD — 



CD 

o 







2 
0.0 



o 



o 



A< 



O 
A 



A 



O 



O 



o 



X 



+ 



X 
o 
X 



X 



X 



A X+ 

O A 
XA 

+ 
+ 

X 

+ 

A< 



A 



A 

+ 



+A 



- 1 -*A + A - 



+ 



"XX 



AA * 
A+ + + + 

+ + 



0.2 



0.4 0.6 
EW (1 1 .3yum 



0.8 



.0 



Fig. 7. — Ratio of total flux in [Nelll] 15.5 /im to total flux in PAH 11.3 fim is compared 
with EW of PAH 11.3 /iin (in /im), using data for starbursts from Tables 2, 4, and 6. 
Triangles are Markarian starbursts from Figure 1, diamonds are GALEX starbursts from 
Figure 2, and c rosses are Sp i tzer s tarbursts from Figure 3. Asterisks are nearby starburst 
galaxies in the iBrandl et al.l (120061 ) sample measured with high resolution spectroscopy in 



Bernard-Salas et al.l ( 120091 ). Small symbols show limits on undetectable [Nelll] for all cate- 



gories. Plot includes only starbursts; the AGN shown with different symbols in Figure 4 are 
not included. [NeIII]/PAH measures hardness of the ionizing radiation, and EW(11.3 /an) 
shows strength of dust continuum relative to strength of PAH. 
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Fig. 8. — Total flux of [Nelll] relative to total flux of PAH 11.3/im compared with luminosity 
vL v (7.7 fim) in erg s _1 (log [i/L„(L Q )] = log [i/L u (erg s -1 )] - 33.59). Triangles are Markarian 
starbursts from Figure 1, diamonds are GALEX starbursts from Figure 2, and crosses are 
Spitzer starbursts from Figure 3. Small symbols show upper limits for [Nelll] in objects 
where this could not be measured. The AGN shown with different symbols in Figure 4 
are not included. There is no trend with luminosity for this ratio among the samples, but 
the values show as in Figure 7 the systematically stronger [Nelll] compared to PAH in the 
UV-selected Markarian and GALEX samples compared to the Spitzer sample. This result 
is also evident in the average spectra of Figure 5. 
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Fig. 9. — Star formation rates from Tables 2, 4, and 6 determined from isL u (7.7 /im) relative 
to SFRs determined from z/Lj,(153 nm), compared with luminosity isL u (7.7 jjxa) in erg s -1 (log 
[z/L v (L Q )] = log [^Lj,(erg s -1 )] - 33.59). Triangles are Markarian starbursts from Figure 1, di- 
amonds are GALEX starbursts from Figure 2, and crosses are Spitzer starbursts from Figure 
3. The AGN shown in Figure 4 are not included . Line is the fit to the full Spitzer-selected 
sample of 287 sources in lSargsyan and Weedmanl (120091 ); sources in the Spitzer sample from 
Figure 3 shown with crosses are the most extreme in SFR(PAH)/SFR(UV). The Markarian 
and GALEX UV-selected samples show similar values of SFR(PAH)/SFR(UV), with median 
log SFR(PAH)/SFR(UV) = 0.8. If the discrepancy between SFR(PAH) and SFR(UV) is 
caused only by the obscuration of the UV, this result for ultraviolet-selected sources indicates 
that 15% of the intrinsic UV luminosity if observed. Much greater obscuration is implied 
for the Spitzer sample, with median log SFR(PAH)/SFR(UV) = 2.4, indicating that only 
0.4% of the intrinsic ultraviolet luminosity emerges. 
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Fig. 10. — Star formation rates from Tables 2, 4, and 6 determined from vL l/ (7.7 /im) relative 
to SFRs determined from z/Lv(153 nm), compared with luminosity as measured by M(UV); 
M(UV) is an AB magnitude determined at rest frame 153 nm as explained in Tables 2, 4, 
and 6. Triangles are Markarian starbursts from Figure 1, diamonds are GALEX starbursts 
from Figure 2, and crosses are Spitzer starbursts from Figure 3. The AGN shown in Figure 4 
are not included. Line shown is the formal linear least squares fit to all of the points, of form 
log[SFR(PAH)/SFR(UV)] = (0.29± 0.04)M(UV)+6.67±0.67. This line has the opposite 
slope of the line in Figure 9, which shows that objects selected using UV selection criteria 
favor sources having less obscuration. Points show the similarity of Markarian and GALEX 
sources in terms of obscuration as measured by SFR(PAH)/SFR(UV), but show that the 
GALEX sample is systematically more luminous. 
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Fig. 11. — Observed flux density of 7.7 PAH feature, f u (7.7 jum), relative to observed 
ultraviolet flux density from GALEX, f„(153 nm), from Tables 1, 3, and 5 compared to 
luminosity i/L„(7.7 /im) in erg s _1 (log [i/L„(L )] = log [^L„(erg s -1 )] - 33.59). Triangles are 
Markarian starbursts from Figure 1, diamonds are GALEX starbursts from Figure 2, and 
crosses are Spitzer starbursts from Figure 3. The AGN shown in Figure 4 are not included. 
This illustrates the empirical result that the smallest value of log[f;,(7.7 //m)/f„(153 nm)] = 
1. We take this as the empirical value which relates PAH luminosity to ultraviolet luminosity 
for a starburst with negligible obscuration of the UV. Assuming no obscuration at 7.7 /im , 
the fraction of UV luminosity intrinsically emitted compared to the luminosity observed, 
UV(intrinsic)/UV(observed), is then given by log[UV(intrinsic)/UV(observed)] = log[fj,(7.7 
//m)/fj,(153 nm)] - 1. 
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Fig. 12. — Observed flux density of 7.7 PAH feature, f v (7.7 jum), relative to observed 
ultraviolet flux density from GALEX, f„(153 nm), from Tables 1, 3, and 5 compared to 
luminosity M(UV). M(UV) is an AB magnitude determined at rest frame 153 nm as explained 
in Tables 2, 4, and 6. Triangles are Markarian starbursts from Figure 1, diamonds are 
GALEX starbursts from Figure 2, and crosses are Spitzer starbursts from Figure 3. The 
AGN shown in Figure 4 are not included. Trend of points is the opposite of trend in Figure 
11, which shows that the most luminous sources in M(UV) are also the least obscured. This 
figure also shows the empirical result that the limiting value of log[f^(7.7 /xm)/f I/ (153 nm)] 
= 1, which is taken empirically to be the value for an unobscured source. 
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Fig. 13. — Ratio of intrinsic ultraviolet luminosity to emerging luminosity compared 
to M(UV) for all samples of starbursts. Triangles are Markarian starbursts from Fig- 
ure 1, diam onds are GALEX starbu r sts fro m Figure 2, and crosses are all Spitzer star- 
bursts from ISargsyan and Weedmanl ( 120091 ). The AGN shown in Figure 4 are not in- 
cluded. Intrinsic(UV)/observed(UV) is determined from relation in section 5.2, that 
log[UV(intrinsic)/UV(observed)] = log[f v (7.7 //m)/f I/ (153 run)] - 1. Lines shown are 
linear least squares fits to the different samples. The fit to the ultraviolet-selected 
Mar karian-l- GALEX sample (lower line) has the form log[UV(intrinsic)/UV(observed)] 
= 0.07(±0.04)M(UV)+2.09±0.69; the fit to the full infrared-selected Spitzer sample 
(central line) is log[UV(intrinsic)/UV(observed)] = 0.17(±0.02)M(UV)+4.55(±0.4); the 
fit to the extreme infrared-selected Spitzer sample in Table 5 (upper line) is log 
[UV(intrinsic)/UV(observed)] = 0.14(±0.03)M(UV)+4.68±0.50. 
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Fig. 14. — Ratio of fluxes Ha/H{3 from SDSS compared to the f u (7.7 / um)/f l/ (153 nm), from 
Tables 1, 3, and 5. Relative line fluxes are measured using AfA taking fA at the peak of the 
line. Triangles are Markarian starbursts from Figure 1, diamonds are GALEX starbursts from 
Figure 2, and crosses are Spitzer starbursts from Figure 3. The AGN shown in Figure 4 are 
not included. Intrinsic ratio in unobscured source should be 2.87, which agrees with lowest 
value on plot and confirms determination in Figures 11 and 12 that log[f„(7.7 /xm)/f„(153 
nm)] = 1 for unobscured source. Line shown is linear least squares fit to all of the points, 
of form f(Ha)/f(H/3) = (2.06±0.21)log[f„(7.7 fan) /f„ (153 nm)]+0.93±0.53. The slope shows 
that there is a relation between reddening in the Balmer lines and the obscuration shown by 
comparing PAH and UV flux densities. 
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Fig. 15.— Obscuration at Ha determined from reddening of Ha/H/3 compared to ob- 
scuration of ultraviolet determined from i u (7.7 //m)/f„(153 nm). Triangles are Markarian 
starbursts from Figure 1, diamonds are GALEX starbursts from Figure 2, and crosses are 
Spitzer starbursts from Figure 3. The AGN shown in Figure 4 are not included. Line is linear 
least squares fit to all of the points, of form log [intrinsic Ha/observed Ha] = 0.32(±0.03)log 
[intrinsic UV/observed UV] + 0.16±0.04. Line shows that much more obscuration is present 
in the ultraviolet than would be indicated by obscuration at Ha. 
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Fig. 16. — Ratio of starburst luminosity to luminosity of stellar continuum shown by ratio 
f^(7.7 yum)/f l/ (1.6 /im) compared with luminosity i/L„(7.7 jum) in erg s _1 (log [i/L„(L Q )] = 
log [z/Lv(erg s -1 )] - 33.59). Triangles are Markarian starbursts from Figure 1, diamonds 
are GALEX starbursts from Figure 2, and crosses are Spitzer starbursts from Figure 3. 
Asterisks are very dusty, luminous starbursts at z ~ 2 discovered by Spitzer with f„(24 
jum) > 0.5 mJy The AGN shown in Figure 4 are not included. Small symbols are sources 
undetected in 2MASS, for which a limit is shown corresponding to H > 16.5. The starburst 
to stellar ratio is similar among all of the samples and not dependent on luminosity. 
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Fig. 17. — Total bolometric luminosity L ir in erg s _1 compared to M(UV). (log [I/j r (L )] 
= log [ Lj r (erg s -1 )] - 33.59.) Triangles are Markarian starbursts from Figure 1, dia- 
monds are GALEX starbursts from Figure 2, and crosses are Spitzer starbursts from Figure 
3. Luminosities Lj r are determined from the empirical transformation log L ir = log[z/L„ 
(7.7/Lim)] + 0.78 (IHouck et al.l 120071 ) using the ^L v (7.7/im) from Tables 2, 4, and 6. Up- 
per line is linear least squares fit to extreme Spitzer sample in Table 5, of form log Lj r = 
-(0.25±0.03)M(UV)+40.79±0.48, and to the GALEX+Markarian ultraviolet-selected sam- 
ples in Tables 1 and 3, of form log L ir = -(0.33±0.04)M(UV)+38.10±0.69. Lines represent 
the extremes of infrared selection and ultraviolet selection for starbursts and indicate that 
ultraviolet-selected starbursts of a given M(UV) underestimate bolometric luminosities by 
about a factor of 30 over all luminosities when compared with the most obscured starbursts. 
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Fig. 18. — Total bolometric luminosity L ir in L Q compared to M(UV). (log [L ir (L )] = 
log [Lj r (erg s -1 )] - 33.59.) Triangles are Markarian starbursts from Figure 1, diamonds 
are GALEX starbursts from Figure 2, and crosses are Spitzer starbursts from Figure 3. 
Luminosities L ir are determ ined for all source s from the empirical transformation log L ir 
= log[i/L„(7.7/im)] + 0.78 fjHouck et al.l 120071 ) using the vL v (7.7/im) from Tables 2, 4, 
and 6. Lines shown are fit separa tely to the full Spitzer-selected sample of 287 starbursts 



m 



Sargsvan and Weedmanl (J2009J), of form log L ir = -(0.23±0.02)M(UV)+6.99±0.41 (up- 



per line), and to the GALEX+Markarian ultraviolet-selected samples, of form log 
-(0.33±0.04)M(UV)+4.52±0.69 (lower line). 
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Fig. 19. — Ratio log L ir /\og[uL u (7 .7 nm)] compared to luminosity vL„ (7.7 urn) in e r g s" 1 . 
(log [z/L„(L )] = log [z/Lj,(erg s™ 1 )] - 33.59.) Asterisks are starbursts from lBrandl et al.l (120061 ) 
used to determine empirical transformation which we use, log Lj r = log[z/L^ (7.7/im)] + 0.78. 
Vertical lines are results for Markarian starbursts in Table 2. Results are uncertain because 
of uncertain aperture corrections when comparing IRS and IRAS flux densities. Upper value 
for each line is value found if total IRAS flux is compared to observed i u (7.7 /mi). Lower 
value for each line is value found if IRAS fluxes are all decreased by a factor corresponding 
to the ratio f^IRS 25 /xm)/f l/ (IRAS 25 fim) explained in Table 2. Small triangles represent 
upper limits for sources without IRAS detections in one or more bands. 
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Fig. 20. — Total bolometric luminosity L ir in L compared to M(UV) for samples discussed 
in present paper and shown in Figu res 17 and 18 (upper t hree lines) compared with summary 
conclusions for LBGs at z = 2.5 in lBouwens et al.l (120091 ) (line with asterisks). Upper line is 
for the extre me Spitzer starbursts in Table 5; second line is for all of the infrared-selected 
starbursts in lSargsyan and Weedmanl (120091 ). and plotted crosses show this full sample; third 
line is the combined Markarian and GALEX ultraviolet-selected samples in Tables 1 and 3. 
Results indicate that the adopted result for LBGs at high redshift (line with asterisks) applies 
only to the lower envelope for all starbursts. Differences between this line and other lines 
show differences that would arise for measures of SFR density. Adopted results for LBGs 
agree with the full infrared-selected sample at M(UV) ~ -21, but differ by a factor > 10 at 
M(UV) = -17. Differences are greater at fainter M(UV) for reasons discussed in section 5.2, 
that the faintest sources in the ultraviolet are those with the greatest obscuration by dust. 



